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We describe fermions in terms of a classical statistical ensemble. The states r of this ensemble 
are characterized by a sequence of values one or zero or a corresponding set of two-level observ- 
ables. Every classical probability distribution can be associated to a quantum state for fermions. 
If the time evolution of the classical probabilities pr amounts to a rotation of the wave function 
qrit) = ±^/pAf), we infer the unitary time evolution of a quantum system of fermions according 
to a Schrodinger equation. We establish how such classical statistical ensembles can be mapped 
to Grassmann functional integrals. Quantum field theories for fermions arise for a suitable time 
evolution of classical probabilities for generalized Ising models. 
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I. INTRODUCTION 

It is widely believed that fermions are a genuine quan- 
tum concept that is not present in classical statistics. In 
contrast, we establish here that fermions arise rather nat- 
urally from a classical statistical ensemble. This is not in 
contradiction with the quantum features of fermions. Ac- 
tually, all the usual quantum features arise together with 
the fermions from suitable classical statistical ensembles. 
We encounter here an example for the emergence of quan- 
tum physics from classical statistics [1]. 

A classical statistical ensemble is characterized by the 
states r and the normalized positive probabilities Pr for ev- 
ery state, Pr > 0, Pr = 1- We consider here a set of two- 
level observables Na which take in every state r the values 
zero or one. We further assume that the states can be com- 
pletely labeled by these values Ua = 0, 1. This means that 
r = [rifj] can be specified by an ordered set of numbers 
or 1 for the B two-level observables, a, /3 = 1 . . . B. For 
finite B there are 2^ states, and we also consider the limit 
B ~> oo. This setting is a very general one, since arbitrary 
observables can be composed from two-level observables. 
For example, the position observable of a particle can be 
defined by dividing space into cells by yes/no questions as 
upper half or lower half, left or right etc., and making such 
a grid as fine as required by a given experiment. Each 
yes/no question can be associated with a two-level observ- 
able. Already at this early stage the possible association 
between the set of numbers [np] and the occupation num- 
bers for a quantum system of fermions is suggestive. 

If a consists of a position label x, for example denoting 
the points of a D-dimensional hypercubic lattice, plus an 
"internal index" 7 = 1.../, a = (a;,7), the states t are 
the possible states of a classical Ising model with / different 
species. Instead of the occupation numbers ria we can also 
use unit spins Sa = 2na — 1 = ±1. In this case every 
probability distribution {pr} denotes a generalized Ising 
model. (A particular family of {pr} corresponds to the 
standard Ising model with next neighbor interactions.) We 
will associate to each {pr} a quantum state for a multi- 
fermion system. In this sense we map classical statistical 
ensembles of a generalized Ising type to quantum fermions. 

A crucial concept for the description of fermions is the 



"classical wave function" {gr}, which consists of the posi- 
tive or negative roots of the probabilities Q 



(1) 

As we have discussed earlier 3 : 0] , the relative signs of the 
real numbers qr are largely fixed by continuity properties, 
such that only few relative signs without physical relevance 
remain free. If all physical quantities are ultimately com- 
putable in terms of the probability distribution {pr}, the 
signs in the classical wave function {qr} carry no indepen- 
dent relevant information. The choice of these signs can be 
associated with a choice of gauge. We can then use {qr} or 
{pr} as equivalent descriptions of the classical statistical 
ensemble. For all observables and correlations built from 
Na the sign s,- = ±1 of (7,-, 



fPr, 



(2) 



obviously plays no role, since expectation values can be 
computed from {pr}- 

One may further introduce a large class of " statistical 
observables" Q. They can be represented as off-diagonal 
operators acting on qr. For those the signs Sr matter, 
but they are fixed in a way such that the statistical ob- 
servables remain computable in terms of {pr}- Again, for 
such observables all necessary information is contained in 
the probability distribution {pr}- In the present paper the 
statistical observables will not play a central role and we 
rather concentrate on the "diagonal" observables built from 
Na which take a fixed value Ar for every classical state r. 
Nevertheless, the use of the classical wave function {qr} is 
very convenient for the description of the classical statisti- 
cal ensemble and constitutes the key for its interpretation 
in terms of fermions. (There are some common features 
between the classical wave function and the Hilbert space 
formulation of classical mechanics by Koopman and von 
Neumann Q. In contrast to this work the classical wave 
function is real, however, and this will be a crucial ingre- 
dient for our approach.) 

In this paper we concentrate on the description of the 
ensemble in terms of the classical wave function. We will 
identify {qr} with a quantum wave function. In the case 
of Ising like systems with a = {x, 7) this wave function 
describes a quantum theory for an arbitrary number of 
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fermions. We further map each wave function {qt} to an el- 
ement of a Grassmann algebra. This allows the implemen- 
tation of the fermionic creation and annihilation operators 
and the discussion of quantum states with a fixed number 
of fermions. We discuss symmetries that can be realized in 
a simple way on the level of fermions, while on the level of 
the classical probability distribution {pr} they are rather 
hidden in the form of non-linear transformations. All this 
discussion is independent of the particular dynamics of the 
system. The first part of our paper therefore constitutes 
a very general map from a discrete classical statistical en- 
semble to a quantum theory of fermions. In this respect it 
is in line with earlier implementations of fermions within 
bosonic systems 3. 

The second part of this paper concerns the time evolu- 
tion of the probability distribution of the classical statisti- 
cal ensemble. For this purpose we need an evolution law 
which describes how {pr} depends on time t. This evolu- 
tion law has to be specified as the basic dynamical law - it 
is not known a priori. There is substantial freedom in the 
choice of such a law - one only has to guarantee that all 
probabilities stay positive at all time, Prit) > 0, and that 
the normalization is preserved, '^T-Prit) = 0. 

We will discuss a specific form of the time evolution of the 
probability distribution as given by the "unitary evolution 
law" 

==^i?,pgp(i-e), (3) 
p 

with rotation matrix R independent of {qt}- (The most 
general evolution of {pr} can be described with R depend- 
ing on {qt}-) The evolution law allows the computation of 
the classical wave function q^- (t) at a time t in terms of the 
wave function {t — e) at a previous time t — e which may 
be infinitesimally close. The form of a rotation preserves 
the length of the real unit vector q — {qr} and therefore 
the normalization of the probability distribution, 

q' = Y.<lr=T.Pr = ^- (4) 
r r 

It is easy to describe periodic changes of {pt} in terms of 
such rotations 2]. 

In the limit e — > we may cast the rotation of q in eq. 
([3]) into the form of a differential equation 

dtq^Kq , K ^ {dtR)R^^ , = -K. (5) 

For any time evolution of q we can compute the time evo- 
lution of the probability distribution {pr} = {q^}- While 
the evolution equation ([5|) for g is a linear differential equa- 
tion, the corresponding differential equation for {pt} can 
be non-linear. We may use the standard formalism of quan- 
tum mechanics by defining K — —iH/h. In our approach 
the value of Planck's constant h is purely an issue of units. 
(The classical limit of quantum mechanics is realized as 
usual if the action is large in units of h.) Using inverse 
time units for H we will set H = 1. The evolution is then 
described by a familiar type of "Schrodinger equation" 

dtq = —iHq. (6) 



For finite B the Hamiltonian H is a purely imaginary 
antisymmetric 2^ x 2^ matrix and therefore hermitean, 
W = H . At this point we have already the formal struc- 
ture of the unitary time evolution of a wave function in 
quantum physics. We will see for detailed examples how 
such classical ensembles describe quantum systems. So far 
the wave function is real. We will explain in our specific ex- 
amples how the complex structure characteristic for quan- 
tum physics arises. We will concentrate here on the de- 
scription of quantum systems for fermions, while bosonic 
objects can be found as composites involving even numbers 
of fermions. 

Starting from these rather simple and very general con- 
siderations we will see how the whole formalism for a quan- 
tum field theory of multi-fermion systems arises from a 
classical statistical ensemble. For the time evolution of the 
system we exploit the isomorphism between wave functions 
{q-r} and elements of a Grassmann algebra. We construct 
a map between a family of wave functions {9r(i)} for ev- 
ery time t on one side, and a Grassmann functional in- 
tegral on the other side. This allows us to use all tools 
of quantum field theory for the description of a unitary 
time evolution of classical probabilities ([3]). In particular, 
we formulate the time evolution for an Ising type classical 
statistical model in two dimensions which is equivalent to 
a two-dimensional Lorentz-invariant quantum field theory 
fermions. 

This paper is organized as follows. In sect. |ll]we intro- 
duce the Grassmann representation of the wave function 
{qr}- Up to irrelevant signs this also constitutes a Grass- 
mann representation of the "classical" probability distri- 
bution {pt}- In sect, mil we briefiy recapitulate the con- 
nection between the Grassmann representation of the wave 
function and the quantum wave function of a multi-fermion 
system in Fock space. Sect. IIVI is devoted to conjugate 
Grassmann variables. We show how arbitrary Grassmann 
operators can be represented as multiplicative operators 
constructed as linear combinations of powers of the Grass- 
mann variables ip and their conjugates t^. In sect. IVl we 
discuss symmetry transformations of the probability distri- 
bution and wave functions. 

The second part of this paper establishes the map be- 
tween Grassmann functional integrals and unitary time 
evolutions of classical statistical ensembles. In sect. IVII 
we construct the wave function {(j'r(i)} for every time t 
by a partial integration of the degrees of freedom of the 
Grassmann functional integral. We first use a very simple 
example corresponding to two-state quantum mechanics. 
In sect. IVIII we indicate the classical probability distribu- 
tion and its time evolution which corresponds to the func- 
tional integral and quantum wave function of sect. IVII The 
opposite direction is followed in sect. I Villi where we start 
from the time evolution of classical probabilities and con- 
struct the associated Grassmann functional integral. We 
generalize our example and discuss the constraints that the 
Grassmann functional integral has to obey in order to per- 
mit a one to one correspondence with a unitary evolution 
of the wave function {gr}- 

While the first example in sects. IVIIIVIIII emplovs a real 
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Grassmann algebra, we present a simple example for a 
complex structure and the resulting complex Grassmann 
algebra in sect. IIXI Sect. |X] addresses the issue of the 
boundary terms in the functional integral which specify 
the specific quantum state. If one of these terms is formu- 
lated as an initial condition for the past, while the other 
concerns the future, we need to relate the future to the 
past. This can be done using translations and refiection 
in time. Finally, in sect. IXII we discuss a two-dimensional 
quantum field theory for free fermions within our setting. 
We derive the time evolution of the classical probability 
distribution {prit)} from a Lorentz-invariant Grassmann 
functional integral. Inversely, for a suitable time evolu- 
tion of the classical probability distribution for an Ising 
type system we construct the quantum wave function for 
a multi-fermion system and its time evolution. This can 
then be mapped to a functional integral, from which the 
observables in a quantum field theory can be computed in 
the usual way. This system admits a complex structure 
which can be associated to the notion of antiparticles. We 
present our conclusions in sect. I XIII 



II. GRASSMANN REPRESENTATION OF 
CLASSICAL PROBABILITIES 



We want to show that classical statistical ensembles with 
probability distribution {pr} following a "unitary time evo- 
lution" ([S]) can describe quantum systems for fermions. In 
particular, we want to recover the standard formalism with 
creation and annihilation operators. For this purpose we 
first map in this section the classical wave function {qt} 
to an element of a Grassmann algebra, and connect the 
two-level observables Na to operators acting within the 
Grassmann algebra. The connection between the Grass- 
mann algebra and fermionic quantum systems is standard 
and will be briefiy recapitulated in the next section. In 
sect. IXII we will construct a classical statistical ensemble 
describing a quantum field theory for fermions. From this 
we can easily obtain quantum field theories for bosons as 
suitable composite fields for fermion bilinears or more com- 
plicated bosonic composite fields. This may include gravity 
and gauge interactions as in spinor gravity . 

The map from the classical wave function to an element 
of a Grassmann algebra has been constructed in detail in 
ref. [2}. We briefly recapitulate in this section the repre- 
sentation of states, probabilities, observables and evolution 
law in terms of a Grassmann algebra. The two-level observ- 
ables Na can be associated with bits taking the values one 
and zero. To each bit we associate a Grassmann variable 
ipajCX — 1...B. Grassmann variables anticommute and 
obey the standard differentiation and integration rules 



(7) 



dip/3 



To every state r = [n^] we associate an element of the 
Grassmann algebra, t ^ g^, which is a product of factors 



tpp according to the following rule: if within the sequence 
[n^j] one has for a given a the value — 0, we take a factor 
tjja in Qt, while na = 1 corresponds to a factor 1. Thus an 
ordered chain (ni,n2, . . . ,713) corresponds to a product of 
factors ipa for each tIq = in the chain. The Grassmann 
elements g^- form a complete basis of the real Grassmann 
algebra G- 

To every wave function {qr} ~ {sr^/Wi-} we associate an 
element of the Grassmann algebra g (z Q, 



(8) 



We also define the conjugate element g of g by 

5 = X] '^■^9r (9) 

r 

where gr obtains from r = [n/3] by taking a factor ipa for 
every = 1, and a factor 1 for = 0, with an appropri- 
ate sign such that grdr = ipiip2 ■ ■ ■'4'b- More formally, we 
may define gr by the property 

j 'Dipgrgp = 5rp, (10) 

where the integral over all Grassmann variables reads 



ViP= diPe--- diP2 dip 



(11) 



The elements g and g associated to the probability distri- 
bution {pt} are normalized according to 

j Vip gg = j Vip^qrgrqr'gr' =^Pt = 1- (12) 

We can realize a classical observable A as an operator A 
acting on g according to 

Agr = Argr , / VlPgAg = ^Pr^r = {A). (13) 

In particular, the occupation number operator A/'q associ- 
ated to Na reads 



Ma = daljja 



d 

dlpa 



Ipa- 



(14) 



Different occupation numbers commute, JVaJ^p = MfjMa- 

The element g — grffr plays the role of a Grassmann- 
valued wave function, similar to a vector which obtains as 
a sum over basis elements gr with coefficients q-r . The time 
evolution of {gT-(i)}, 



9- 



(15) 



is transfered to the time evolution of .g(<), 

g{t) = U{t, to)g{to) , m = ^^{t, toMto)- (16) 
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Here the Grassmann operator U plays a role similar to the 
unitary evolution operator in quantum mechanics, with U 
the inverse of U, 



(17) 



(For more details see ref. 0-) In the Grassmann formula- 
tion the Hamilton operator T-L is defined as 



(18) 



and the time evolution ()16|) can be written as a Schrodinger 
equation 



dtg{t) = -in{t)g{t). 



(19) 



For a real Grassmann algebra the operator /C — —iH is 
real such that the appearance of i in eq. (jl9p serves only 
for the formal analogy to quantum mechanics. 

If a suitable complex structure exists we can map the ele- 
ments of a real Grassmann algebra to elements of a complex 
Grassmann algebra with 2^~^ basis elements. For exam- 
ple, it may be built from F = B ~ 1 Grassmann variables. 
Such mappings will be discussed in our explicit examples. 
They correspond to maps from a real vector {qt} to a com- 
plex vector {cr} with half the number of components. The 



probabilities read now pj- — 



Replacing in eq. © 



Qr — >■ c* the normalization (jl2p remains valid for the com- 
plex Grassmann algebra as well. The time evolution (jisp 
or (|16p will then be mapped to the unitary time evolution 
of a complex wave function {cr}, with Rrp replaced by a 
unitary matrix Urp- Within the description by a complex 
Grassmann algebra H needs no longer to be purely imagi- 
nary. 



III. FERMIONS 

In this section we recall how each element of the Grass- 
mann algebra can be associated with the quantum wave 
function for a multi-fermion system. This material is rather 
standard and displayed here only for completeness and for 
setting the notation for later sections. We will work here 
with a complex Grassmann algebra where the elements can 
be labeled by the components Cr of a complex vector. (The 
case of a real wave function and real Grassmann algebra 
is recovered by using real Cr ~ Qt-) The multi-fermion 
system can be described in the occupation number basis. 
For only one species of fermions the Hilbert space can be 
built from two states, the vacuum |0) (empty state) and 
the one fermion state |1) (occupied state). The occupied 
state obtains by applying the creation operator a' to the 
vacuum 



|l)^at|o), 



(20) 



where the annihilation operator a acting on the occupied 
state yields the vacuum 



a|l) = |0). 



(21) 



Using further the relations 

a|0) = a^|l) = 



(22) 



implies — (a^)^ = 0. The occupation number operator 

M = a)a acts as 



i^a|l) ^ |1), a^a\{)) = 0. 



(23) 



Thus |1) and |0) are eigenstates of M with eigenvalues one 
and zero. Furthermore, with 

aa^\l) = 0, aa^\0) = |0), (24) 

we infer the commutation relation 

{a\a]^l. (25) 

The wave function for a pure quantum state for one species 
of fermions can be written as 



VJ = ci|l)+co|0), |cip + |co|2 = l. 



(26) 



Expressed in terms of a two-component complex vector 
(ci,co) one has the explicit representation 



0, 0\ t / 0' 1 

°= ' 1, ' " = 0, r 



(27) 



The operators and states for one fermion species can 
be identified with the Grassmann operators and elements 
of the Grassmann algebra built from a single Grassmann 
variable ip. We associate 



dip ' dip ' 



(28) 



where the Grassmann operators act on states represented 
by elements of the Grassmann algebra. The basis elements 
are associated as 



|o)=v. , |i)=i 



such that 



V = ci|l) 



Co 



g = cii> + cq. 



(29) 



(30) 



In other words, g is specified by the two complex num- 
bers {ct-} = {ci, Co} which characterize the quantum wave 
function (|26l) . 

The description of the quantum wave function for 
fermions by an element of the Grassmann algebra is eas- 
ily extended to several fermion species that we label by 
a — 1 ... F . We now use a Grassmann algebra based on F 
Grassmann variables i>a. The basis states of the Hilbert 
space (Fock space) for the multi-fermion system obtain by 
applying various powers of creation operators ajj for dif- 
ferent fermions species to the vacuum. For the example 
of two species {F — 2) we have four basis states, the vac- 
uum |0, 0), the two states |1, 0) and |0, 1) with one fermion. 
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either of type 1 or type 2, and the totally occupied state 
|1,1). They obey 

|1,0) = al|0,0) , |0,1)=4|0,0), 
|1,1> = -a\al\0,0) = 4a||0,0) 

= -a\\0,l)=al\l,0), (31) 

and 

|1,0) = a2|l,l) , |0,1) = -ai|l,l) 
|0,0) = aia2|l,l) =-a2ai|l,l) 

= a2|0,l) =ai|l,0). (32) 

Our sign conventions for the basis states are compatible 
with the anti-commutation relations 

{tta, a^} = , {al,al} = , {aj,, a/j} = S^is- (33) 

A pure state quantum wave function reads 

V'^^Crifr, (34) 
r 

with 

VPi = |1, 1) , (^2 = |0, 1) , |1, 0) , (^4 - |0, 0). (35) 

The map to the Grassmann algebra employs the identifi- 
cation 



d 



d 



. - , a„ = , AA„ = Q^J., (36) 



and the association 



|0,0) SVjlV^2 = g4 , |1,0) = V^2 =33, 

|0,1) ^ -Vi=g2 , |1,1) = 1=51, (37) 



with 



f = 9 = > , Cr9 



(38) 



The relations dM]), §6}, §8^ hold for an arbitrary 

number of fermion species F. Our conventions are such 
that states which obtain from the vacuum by applying an 
ordered sequence of creation operators aj^ , with larger a to 
the left, have a plus sign. This concludes the one to one 
correspondence between multi-fermion wave functions and 
elements of a complex Grassmann algebra. 



IV. CONJUGATE GRASSMANN VARIABLES 

The representation of Grassmann operators in terms of 
■0 and d/d'ijj is somewhat cumbersome, since these two ob- 
jects do not anticommute. We therefore discuss in this 
section another representation of the classical probability 
distribution, wave function and observables where d/ dtp is 



effectively replaced by a new Grassmann variable ip, obey- 
ing 

{V'q,'0/3} = {V'a, = {'0a,'0/i} = 0. (39) 

1. Conjugate variables 

This goal will be achieved by a type of Fourier transforma- 
tion of g, 



with 



Vip = I dipi I dip2... / dtp 



(40) 



(41) 



We will systematically replace g in sect. HI] by the associ- 
ated element g{ip). 

The map 5 — > ^ in eq. (j40[) is invertible and we next 
construct the inverse map g ^ g. Basic identities needed 
for this purpose can be found, for example, in ref. Q. For 
general Grassmann variables obeying eq. (1391) we note the 
identities (no sum over /3 here) 

exp(-0^V'/3) = 1 + '0^V'/3, 
dippd'ipfie-yiYii'ipfiipp) = 1, 
d 



dipp 



exp 



i^^p-ipp) = --ipfse-xpi'ipfs-ipp)- (42) 



For a given a, /3 one has a (5-function like object 

^{'Ipa - Ipli) = V'q - V"/? = -<5(?/'^ - -4^0), 

dlpaS{i)a - tpi3)h{lpa) = Hlp/^), (43) 



where h is an arbitrary element of the Grassmann algebra 
which may also involve Grassmann variables different from 
tpa- We can write 



S{ipa - -0/?) = / dipp exp {ip/^itpa - Ipf))} (44) 
dtpp exp{-tpf3tpi3) exp{tpaipp), 



and 



Si^pa-^p) = j di}peii-p{'lpj3{%pa- 

= -6{tPp-^^). (45) 
These identities are easily generalized 

/ T>ip'D%p exp('y^ tppipf^) — 1 

exp(^ -ippipn) = --ipa exp(^ 0/30/3), 



dtpa 



- ^) = (-1)^ j PV^exp { ^(0;, - 0^)0^}, 



vip's{tjy - tp)h{ip') = hiip). 



(46) 
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The identities ([46]) can be used in order to verify that 
the inverse transformation g ^ g obeys 



(47) 



Since eqs. ([5]), ^ (and similar for a complex Grassmann 
algebra) define a map g ^ g, we can combine this with eq. 
(|47l) for the definition of a map g g. In fact, the element 
g{ip) can be easily obtained from g{il') by the following 
steps: (i) each factor ipa in g is replaced by ipa in g, (ii) 
the factors ipa are totally reordered in g, (iii) for a complex 
Grassmann algebra one takes the complex conjugate of the 
coefficients Ct . For example, the element associated to g = 
is 3 = c*ip2'^i- The normalization of the conjugate 
basis elements g-r obeys 



/ I?V;exp(^V;^iA/3)ffp(V')5r(V') = 6pr- (48) 



We also note the useful identities 

X!5r(?A)gr(V') = exp(^ V/3'/'/3), 



X!5r(V')5r(V') = exp(^V);3V/3)- 



(49) 



/3 



In terms of the conjugate basis elements we can write 

.9 = X! '^^5^ ' 5 = X! '^^^r, (50) 

r r 

with Ct- = c* = (7t- for a real Grassmann algebra. 

2. Conjugate Grassmann operators 

We next define for an arbitrary Grassmann operator J- the 
normal ordered form J"Ar. It consists of sums of products 
of factors a + htpp + cd/dtpp + d{d / dip p)tp p , i.e. where 
d/dipp always precedes a possible factor ipp. We wiU re- 
place FN{d/dtpi3,ipp) by FN{tpp,^p) according to 

V^gi^)J^N (^''^)5(V') 

2?V2?^exp(^V^;3V^)g(V')i^Ar(V',V')ff(V')- (51) 

Here F^lipjip) obtains from J'n (^■^i'>P^ by replacing 

8/ dip — >■ ?/;. The proof of the identity ([5T|) can be sketched 
as follows. By partial integration we find for an arbitrary 
element / of the Grassmann algebra 

2?#(V')^/(V') 

2?V2?Viexp(^Vi^V/^)g(iA)^Aa/(V')- (52) 

13 



The identity ([52]) also holds if / takes the form f{ip) = 
(af^'^i^) '^ith s-n arbitrary Grassmann opera- 

tor expressed in terms oi ipjs and derivatives 8/ dip p. Fur- 
thermore, we generalize eq. (15^ to 

VipVij] exp(^ ipf3'ijjp)h{ip, ip)j—f{ip) 



2?V2?Viexp(^Vi0V/3)/i(V'>)V'a/(^) (53) 

/3 



for all elements h obeying 8h/8ipa — 0. Using the anti- 
commutation relation 



8 

8lpa 



(54) 



we can bring any arbitrary Grassmann operator T into its 
normalized form Tj^j-. From the identities (|52|) . ([53)) for 
these ordered operators we finally obtain the relation (|51l) . 

In particular, the number operators (jl4p read in this 
representation 



Na = Ipai'o 



(55) 



mi 



while polynomials of different Afai,J^a2 6tc. are repre- 
sented by corresponding polynomials of ipaiipai, "002 "002 
etc. (Note, however, that J\f^ is not represented by 
(ipaipaf — since the normal ordered form reads 
)n — AAq.) The annihilation and creation operators 
are now represented by ip and ip. 
In conclusion, we have arrived at a representation of the 
classical probability distribution and observables where the 
probability distribution {pr} is represented by the conju- 
gate pair of Grassmann elements g{ip) and g{ip), and the 
observables are polynomials in ip and ip. The elements g, g 
are normalized according to 



D4,giiP)gi4,) = 1, 
j DiP^ jv^j j V4,exp{Y.iPpiPfi). (56) 

They encode the information about the classical ensemble 
{Pt} according to 



Cr = / Dlpgr{'ljj)g{'ljj) , Pr = \Cr\" 



(57) 



For a time dependent probability distribution {prit)} the 
coefficients Cr{t), and therefore g{t) — '^^Cr(t)gr., depend 
on t. 

A classical observable A is represented by an element 
of the Grassmann algebra A{ip, ip) that does not contain 
derivatives anymore. For example, we may investigate the 
(equal time) correlation function of the classical statistical 
ensemble. 
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with Na^^T = the vahie of the observable Na^ for the 
state T. In terms of the Grassmann algebra this correlation 
function is given by 

^i?.L„ (t) = I ^#(^; t)Nc., ■ • . N^^Mi'-, t), (59) 

with Na ~ i^ai^ai providcd that all aj are different from 
each other. 



V. SYMMETRIES 

At a given time t a classical statistical ensemble is char- 
acterized by a probability distribution {pr} which corre- 
sponds to a map t ^ pr- Consider an invertible transfor- 
mation between the states, r' = S{t) , r = S~^{t'). This 
induces a transformation in the space of statistical ensem- 
bles, {pr} {p'r} = S{{pr}), given by p'^, = Pr{r') = 

Pg-i(^^,y The transformation 5* is a symmetry if the prob- 
ability distribution remains uncharged, p'^ = Pt- In this 
case the transformed observable A' — S{A), which takes 
for every state r the value A'^ — A^i(^^'f ~ ^s(r)i has the 
same expectation value as A, 

{A') = Y,KPr=Y.^-'i-)P- 

T T 

r' r' 

In the presence of an evolution law the notion of symme- 
try can be extended to relate probability distributions at 
different times, 

K(ti)} = {Prito)} ^ {A'ih)) = {A{to)). (61) 

We will be interested in symmetries where S{t) is re- 
alized by simple transformations of the occupation num- 
bers TT-Q. The transformation of "charge conjugation" maps 
C{na) = fia = ^ — n^, inducing C(r) = f where f obtains 
from T by flipping all bits. A time reversal combines this 
transformation with a reflection in time around some point 
io, Pf{to + i) = Prito — t). Symmetry under charge conju- 
gation C or time reversal T implies {C{A){t)) = {A{t)) = 
{A{t)) or {T{A){to + t)) = {A{to + 1)) = {A{to - t)) with 
A.^ ^ A.^ . 

In the Grassmann formalism the map r S{t) is re- 
flected by a map of the basis elements g-^ ^ g'^ — igr' — 
iffS(r) = Sr'gr'- The iuduccd map for an arbitrary Grass- 
mann element g ^ g' = X^t 9t"5t '^^^ equivalently be ex- 
pressed by a transformation of the wave function {qr} at 
fixed gr, 

S{g) = 9 =YqTg'r = Yq'r9T, (62) 

r r 

with q'^i = Sr'qT{T')- We can represent the transformation 
5 by a matrix multiplication 

9t=Y1 ^'^p'^p 
p 



where S has only elements 1,-1 and 0, with precisely one 
value 1 or —1 in each row and column, SrpSrp' — \Srp\Sppi. 
With [StpY — \Stp\ tbis guarantees 

p'r^Y. \^^P\PP- (64) 

p 

For a fixed wave function the basis elements transform as 

p p 

A symmetry is realized if {qr} is an eigenvector of S with 
eigenvalue one. This has to hold for a suitable choice of 
signs Si-'- 

The transformations of r which are induced by transfor- 
mations of the occupation numbers Ua may be expressed as 
transformations of the Grassmann variables Consider 
first a transformation of the type 

Sii^c) = V-a = ^ W^p^Pp , W^W = 1. (66) 

This defines a transformation S^gr) — g'r of the basis el- 
ements by replacing in g^- each factor ^/iq, by -0^. The re- 
sulting matrix Srp in eq. (j65l) is orthogonal, which guaran- 
tees that the normalization of the probability distribution 
is preserved. For the particular case where all matrix ele- 
ments equal 1, —1 or 0, with only one element 1 or —1 
in each row or column, this property is transfcrcd to the 
matrix St p. We recover the map of classical states t ^ t' , 
as induced by n'^ = Y.p Wap\ni3. 

We emphasize, however, that every transformation (|66p 
with orthogonal W induces a transformation 5(.gr) and 
corresponding transformations ^({gr}) and 5({pt-}). In 
particular, if >5'({p7-}) = {Pt} we find symmetry transfor- 
mations of the probability distribution that go beyond the 
ones realized by a mapping of states r — ^ r'. For exam- 
ple, the system may realize a symmetry originating from 
continuous rotations of the vector {"iAa}- No correspond- 
ing simple construction as a "rotation among states" is 
available on the level of the occupation numbers tIq,. The 
Grassmann formalism can be a powerful tool for realizing 
simple symmetries of a classical statistical ensemble. More 
generally, the group of all orthogonal matrices Srp in eq. 
(1631) . which leave a given wave function {qr} invariant, may 
be considered as symmetries. On the level of {pr} this may 
correspond to non-linear transformations. 

Furthermore, we may investigate transformations that 
map Grassmann elements to conjugate ones. An example 
is the charge conjugation which transforms 

Clpa = Ipa, (67) 

where ■0q is the conjugate Grassmann element of ipa, with 

J Vlpt/Jalpis = SafI- (68) 

The basis elements are mapped by charge conjugation as 
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A charge conjugation invariant state obeys = qr for an 
appropriate choice of the sign of g'^ . It reaUzes a symmetry 
of the classical probability distribution {p%} = {(g^)^} = 
{q^} = {pt}- For a classical observable A constructed from 
the two-level observables (e.g. equal time classical correla- 
tion functions) the charge conjugate observable A'^ — C{A) 
obtains by replacing each factor Na by 1 — ■ The associ- 
ated Grassmann operator A'^ obtains from A by replacing 
ipa d /d^a- For a charge conjugation invariant state one 
finds {A") = (A). 

The concept of symmetries is extended to the formalism 
with conjugate Grassmann variables in a straightforward 
way. Transformations of the type ([SS)) can be associated 
to corresponding transformations of the wave function (|63p 
and therefore to the transformation of the conjugate Grass- 
mann elements 

g'^J2lr9r^J2lr9r- (70) 

r r 

In the second equation we transform the conjugate basis 
elements g-r, keeping the wave function fixed. This general 
transformation rule for g^ can be extended to arbitrary 
transformations (|63)) . including, for example, a charge con- 
jugation of the type ()67p . The conjugate Grassmann vari- 
ables ipa are particular conjugate basis elements and their 
transformation properties are therefore fixed. Transforma- 
tions of the type ((66|) transform basis elements containing 
only one factor -00 into each other 

(71) 

With 5'"'^' — we find the transformation 

P 

As a result, the factor exp(^^ t/'^V'^) appearing in the 
transformations (|40|) . (j47|) or the measure (|56p is invari- 
ant under such transformations. 

VI. QUANTUM WAVE FUNCTION FROM 
GRASSMANN FUNCTIONAL INTEGRAL 

It is well-known that quantum systems of fermions can 
be described by Grassmann functional integrals. In the 
opposite direction, a Grassmann functional integral obey- 
ing certain conditions naturally leads to a wave function 
{qrit)} obeying a unitarity time evolution, and the asso- 
ciated classical probability distribution {prit)}- In this 
section we present a simple example how the unitary time 
evolution follows from a suitable Grassmann functional in- 
tegral. In the following two sections we take the opposite 
direction, starting with a given time evolution of a proba- 
bility density {^^-(i)}, and inferring the associated Grass- 
mann functional integral. We will present further examples 
in sects. IIXI and IxH 



1. Action and functional integral 

For a formulation of a Grassmann functional integral we 
consider a discrete chain of different times t„, with — 
t„ = e and ti„ < tn < tf. In the end we take the limit 
e for fixed Un and tf. For each i„ we employ a set of 
B Grassmann variables ipaitn) and the associated ipaitn)- 
The action S is an element of the real Grassmann algebra 
constructed from ipaitn) and ipa(tn)- As a simple example 
with i? = 1 we discuss the action 

5 = ^L(0 (73) 

t' 

with 

m = Ht')W + e)- V'(i')) + eu;{i^{t') - tPit' + e)). 

(74) 

The sum over t' extends over all i„ within the interval tin < 
tn < tf, where for t' = tf we define L{tf) = —il}{tf)i/j{tf). 
Taking the continuum limit e — > the action takes the form 

j dt[i{t)dtm + ^{iAt)-Ht))]- (75) 

We may define the hermitean conjugate of a Grassmann 
element G[ip, V'] by exchanging ip ip, a, total transposi- 
tion of the order of all Grassmann variables, and a complex 

conjugation of all coefficients, [ip, V^] = ^G* [tj), ij)]^ . 

With this definition we observe that S is antihermitean, 
= —S. This corresponds to a hermitean Minkowski 
action 

SM=iS , e"^ ^e*-^", (76) 

as usual for path integrals in quantum mechanics. 
The partition function is defined as 

Z^J PV(t')2?V'(i')5(V'(i/))T{e-^['^-'^l}ff(VXt.„)). 

(77) 

Here we have introduced a time ordering T which puts fac- 
tors e"^(*") to the left of factors e"^'*") if ta > tb- This 
specification is necessary for our example because L(t') 
contains a term with an odd number of Grassmann vari- 
ables such that factors e"^^*"^ and e"^'^*''-' do not commute. 
If S contains only terms with an even number of Grass- 
mann variables there is no need for such ordering. The 
functional integration extends over all Grassmann variables 

lP{t'), lP{t'), t,n<t' <tf 

J Vip{t')Vip{t') = Y[J d^{t')d^{t'). (78) 

The factor g(ijj{tin)) is a "boundary term" which will corre- 
spond to the initial condition for the wave function qritm)- 
We will see that the functional integral defines a wave func- 
tion gr(i) only ii g(t(;{tf)'j is properly related to g(ip{tin))- 
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This will restrict the functional integrals which describe a 
wave function qr(t). 

We can associate to S the normalized action 

5jv[i/^(i'),V'(i')] =I]^(^')+ln^ = 5 + lnZ, (79) 
t' 

and define G in terms of Sn 

G[^,Vi] =5(V'(i/))T{e-^~[^'*l}g(^(t„)). (80) 
This Grassmann element is normalized according to 

J Vil}{t')V^{t')G[il},i;] = 1. (81) 

We will see below that we can achieve Z = 1 by a proper 
normalization of gin — g{'4'{tin)) , such that no need for the 
shift ([7^ arises. 

In order to have a well defined Grassmann functional in- 
tegral we use a finite number of time steps Un < t' < tf. 
Without loss of generality, we choose here tin — —tf- 
The functional integral will then be defined as the limit 
limit f — > cx)) lim(e — 0), or we may consider e — 
for some finite tf. The Grassmann element G as well 
as Z and other quantities depend on the boundary terms 
gin = 9 {4' {tin)) and cjf = g{^[tf)) and we will discuss the 
appropriate relation between gf and gin below. The action 
(|73p is invariant under a time reflection, accompanied by 
a simultaneous exchange of ■0 and -0, a change of sign of 
tjj, 'ijj, and a total reordering of all Grassmann variables 

(5h^(-t'), -^{-t')]f = Sm'),ij{t')]. (82) 
2. Wave function 

The expectation value of observables at a given time t in- 
volve only the probability distribution at time t. In other 
words, it involves only the information encoded in the 
Grassmann element g{t). In contrast, the functional inte- 
gral ([5T|) involves the Grassmann element G. This belongs 
to a much larger Grassmann algebra which is built from the 
variables tp{t') and V'(t') for all t' . We therefore construct a 
map G g{t) which focuses on the necessary information. 

We start with the functional integral expression for the 
occupation number N(t) 

{N{t))=Z~' I V^{t')V^{t')gsN{t) 

xf{exp(-^[0(i'),^(i')])}5™, (83) 

with 

N(t) ^ mm- (84) 

We want to show that this equals the expectation value of 
the two- level observable N{t) in terms of a wave function 
Qrit). For this purpose we decompose 

s = s< + s>- mm, (85) 



with 

5< = ^L(0, 
t'<t 

5> = ^L(O + ^(t)V'(0- (86) 
t'>t 

This decomposition is chosen such that S'< depends only 
on ijjit' < t),il:{t' < t) and tp{t), while S'> depends on 
ip{t' > t),'ip{t' > t) and tp(t). We will show that qr{t) or 
the associated Grassmann element g{t) obtains by a suit- 
able functional integral over Grassmann variables at t' < t, 
which involves S'<. 

Using the decomposition (1551) we can write 

{N{t)) - J dmdm)N{t)e^^'^^^'^{m))9{m)), m 

with 

g{m) = Z^' J T^^it' < tWit' < t)f{e-^<}g,n, 

g{m) = J Vm > t)Vm > t)9fT{e-^>}, (88) 

and Z^Z^ = Z. This procedure integrates out the future 
and the past and we are left in eq. (|87|) with a Grassmann 
integral involving only ijj{t) and ipit) a-t a given time t. For 
a given gin we choose gf such that g{'tp{t)) according to eq. 

is the conjugate element to g(ip{t)y Then eq. (|57)) 
coincides with eq. ((59)) . We will discuss this issue in more 
detail below and in sect. I Villi 

With this choice of g/ wc have indeed expressed {N{t)) 
in terms of the wave function g{t) 

{N{t)) = J DtPg{t)N{t)g{t). (89) 

Both g{t) and g{t) are explicitly constructed from the func- 
tional integral by eq. (|88|) . defining the maps G — g{t) and 
G — > g{t). Our procedure easily generalizes to a more gen- 
eral form of the action and to more than one species ijja (t) , 
provided the decomposition (|85|) . (|86)) is possible. If the 
action is normalized and contains only terms with an even 
number of Grassmann variables this expression simplifies 

g{t) = [ V^Vm'^^9^n, 
Jt'<t 

m = [ (90) 
Jt'>t 

with an obvious meaning of the functional integrals. 
3. Time evolution 

Let us next investigate the time evolution oi g[t). From the 
definition (|88|) we can obtain g{t + e) by adding an integral 
over the variables il}{t) and i^it), 

dt9{t) = ^{git + e)-git)) 

= \{j dmdme^v{-m]9{m) 

-g{i,{t + e))]. (91) 
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Here we consider g{t) = a{t) + b{t)ip for a fixed Grassmann 
variable that we take a.s ip = ip{t + e). We generalize L{t) 
in eq. dTQ]) to 



L{t) = m{^{t + e) - m) + ^^H{t), 



(92) 



where H{t) depends on ip{t + e) and ■il^{t). For a general 
Hamiltonian H[ip{t),il;{t + e)] the Grassmann algebra is 
extended to a complex Grassmann algebra. For our partic- 
ular example, however, H is purely imaginary, such that L 
remains an element of a real Grassmann algebra. Inserting 
the general expression (|92l) into eq. (I?!]) one obtains 

dtg{t) = -^[jd^t)di;{t)eM-i^H)gS)~g{i^{t + e))] 

^ -I j dmdmH[mMt + e)\9s{t), 

gAt) = exp[m{^P{t) - ^{t + e))}g{i^{t)). (93) 

This equation is exact for H linear in ipjip and gets correc- 
tions which vanish for e — > in the more general case. 

In order to show that eq. ((93)) describes a unitary evolu- 
tion law we next eliminate the conjugate Grassmann vari- 
able T^, reversing the steps in the construction of sect. IIVI 
Using the identity 



diP{t)d'4j{t)ij{t)gs{t) 



d 



d^{t + e) 



d'ijj{t)d'4j{t)g,{t) 



we get the evolution equation 



dtg{t) = -in 



git) 



(94) 



(95) 



and identify with the Grassmann-Hamilton operator. It 
obtains from H by replacing ip — ^ d/d^p, where the ordering 
puts d/dip to the left of For our example one has 



n 



d_ 



(96) 



Eq. (|95|) realizes the unitary time evolution (fT9|) and can 
therefore be associated to a unitary evolution law for clas- 
sical probabilities. We will discuss in sect. IVIII a unitary 
time evolution of classical probabilities which realizes the 
Schrodinger equation 



d,g{t)=^[i^~^]g{t). 



(97) 



Again, this construction is easily generalized. First one 
extracts from the action the Hamiltonian from eq. (I92p . 
with a generalization of the first term by a summation over 
species. Once H is ordered such that all V'a-factors are to 
the left of the ^/^Q-factors one replaces ipa — >• g^— in order 
to obtain Ti. from H . 



What remains to be shown is that the time evolution 
of g{t) and g{t) as defined according to eq. ([55)1 or ([gO]) . 
is compatible with g{t) being the conjugate of g{t) for all 
t. In other words, we want to establish that for g{t) — 
Er'^rit)gr{ip{t)) One obtains g^t) = Er ^WdriHt)) for 
all t, provided this relation holds for some time to- We first 
note the identity 

Mzit) = f Drg{t)g{t) = / I?V(t)2?^(t)5/e-^5m = 1. 



(98) 

If we define g{t) by eq. ([55)1 the time evolution preserves 
Afz independently of g{t) being conjugate to g(t) or not. 
If .^(^in) is normalized with Afz (tin) = 1 we are guaranteed 
that Z = 1 without the need to normalize the action by 
the shift ([7^ . This corresponds to the observation that we 
can obtain Z = 1 by an appropriate rescaling of gin and gf 
in eq. ^ or 

Furthermore, if g{t) is conjugate to g(t) for all t the norm 



(99) 



and Afz (t) coincide. In this case we can guarantee Z = 1 by 
imposing Afp(tin) = 1. The condition Afp{tin) does not in- 
volve any functional integral and can be formulated purely 
as a normalization condition for the coefficients Cr{tin). 
We will show below that eq. ([M]) implies 

dtgit) = -tUit) , dtm = iH'^m, (100) 

where is related to H by the property that for arbitrary 
Grassmann elements g[ip] and f[ip] one has the relation 



DijH^gf = / Di;gnf 



(101) 



For g{t) defined by eq. ([55)1 (and not necessarily conjugate 
to g{t)) we use the expansions 



g 



^r{t)gr , g^Yl ^■rit)9T- 



(102) 



Then eq. (jlOOl) implies for the time evolution of the coef- 
ficients 



dtCr =Y KrpCp , 5tC^ = - ^ CpKp, 

p p 
This follows from the relation 

dtCr = / Dtpdtggr = i / D^pii^ ggr 



(103) 



i j DrgUgr = -J D^gJ^gpKpr- (104) 



We note that eq. (|102p indeed preserves JVz — J^t'^t^t- 
In order to establish eq. (jlOOp we employ 

dtm = --^{git-e)-git)) =dtg{i'(.t-e)) 

dmdmg{m) (los) 

exp { - - e))m}H{t - e). 
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For an arbitrary element f[ilj{t — e)) this yields indeed 



I f D^j{t - e)g{ij{t - e))Hf{iit - e)). (106) 



The condition that g{t) remains conjugate to g{t) for all 
t amounts to c*(t) = Cr{t). From dtc^ ^ J2p^TpC*p we 
infer the condition 



K* — — K 



(107) 



In terms of the Hamihonian H = iK this requires her- 
miticity of H, 



H"^ = H. 



(108) 



Only for an hermitean Hamiltonian the Grassmann func- 
tional integral can describe the time evolution of classical 
probabilities with constant normalization '^T-Prit) = 1- 
For a real Grassmann algebra ii' is a real matrix and must 
be antisymmetric. Eq. (|107l) is obeyed for our example. 
We observe that the evolutions (jlOOp or (|103|) are invert- 
ible if exp {K{t — to)} is an invertible matrix, which holds, 
in particular, for K = —iH, W = H. 

4. Solution of evolution equation 

The formal solution of the evolution equation (^5]) is 
given by the central identity (|88| or ([90)) . provided g/ — 
g{tf). Here the initial condition for the differential equa- 
tion (I102p can be specified at tin by gin = giUn)- Alter- 
natively, we can specify a particular solution by indicating 
g(to) for some arbitrary time Iq. Then gin = g{tin) and 
g{tf), and therefore also gf — g(tf), can be computed from 
the evolution equation. (For gf and gin not related to the 
solution of the evolution equation the functional integral 
([77| yields a transition matrix instead oi Z.) 

The identity ([55)) or (pi)) can be interpreted in different 
useful ways. We may introduce the Grassmann element 



G< =Z-ir{e-^<}5, 



(109) 



which depends on the Grassmann variables '(/'(^' < t)i 
^jj{t' < t) and tpit). We can then interprete g{ip{t)) as 
an integration over "unobservable variables" 



(V'(O) = / 'Di^it' < t)VtP{t' < t)G. 



(110) 



This underlines the character of the local (in time) wave 
function and probability as a subsystem where degrees of 
freedom of a "total system" are integrated out [Tj]. The 
total system may be associated with a probability distri- 
bution for all times, where "probabilistic time" appears 
only as a particular ordering structure . "Unobservable" 
means in this context that observations are performed only 
with local observables. A similar identity holds for g{ijj{t)^ , 

G^=Z-'gff{e-^>} 

g{i;{t)) = / X>V(i' > t}Vi;{t' > t)G>, (111) 



where G> depends on V'(i' > t),'ij;{t' > t) and ii{t). 

Using eq. (186)) we can decompose Z<T'{e~'^<} into a 
product of factors 



^-lg-L(*„) 



(112) 



where the product is over all n with tin l£ tn < t and 
J^j^z^^ = Z^^ . Each individual factor s„ depends on 

i^{tn), i^{tn) and ^{tn+i) = ^{tn + e)- We define formally 
the Grassmann element 



U{tn+l,tn) = / d'lp{tn)e 



(113) 



which depends on ipitn) and ipltn+i)- The operator U 
introduces a new multiplication law which includes the 
Grassmann integral J dip{tn) 



U{tn+l,tn) O g[^{tn)) ^ J dljj{tn)u{tn+l,tn)g {ll^{tn)) 

di,{tn)di^{tn)[e-'-g{'^{tn))]. (114) 



such that [/(t„+i,t„)o(/((?/)(t„)) depends on V'(in-i-i)- Thus 
U can be viewed also as an operator that maps local wave 
functions at t„ onto local wave functions at tn+i- This 
allows us to write g(^ip{t)^ in the form 

git) = g{m) = Y[{U{tn+l,tn)o)g{^j{t,n)), (115) 



where the operator product is defined in analogy to eq. 
(|114l) . ordered such that larger time arguments are on the 
left. 



VII. CLASSICAL PROBABILITIES FOR 
TWO-STATE QUANTUM MECHANICS 

In this section we specify of an explicit unitary evolution 
law for classical probabilities that realizes the evolution 
equation ([^ . It reads 



p^{t) = cos^(cji)pi(0) + sm^{ujt)pa{0) 
— 2scos(wt) sin(cji) ^ pi {0)po (0) , 



(116) 



where we set to = and s — ±1. We note < pi{t) < 1 
and define po{t) = 1 — pi(t), such that the probabilities 
for the two states 0, 1 are normalized, pa{t) + pi{t) = 1. 
If we interprete pi(t) as the probability for the occupied 
state, and pa{t) as the one for the empty state, the mean 
occupation number obeys 



{N{t))^p,{t). 



(117) 



The rotation matrix R in eq. ([3]) can be written as a real 
unitary matrix U 



R 



cosuit, —sinujt 
sinujt, coswi 



U = eiip{-iujT2t). (118) 
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This corresponds to two-state quantum mechanics with 
a hermitean Hamiltonian H = and a real two- 

component initial wave function (f{Q) = (91(0), go (0)) = 
{\/'plW)■,s^/pQ{fi)). Using the standard quantum formal- 
ism for a two-component wave function one has 



-iHip{t) , H = LUT2. (119) 



The wave function ip(t) = U{t,0)ip(0) remains real in the 
course of the evolution. Evaluating the expectation value 
of the operator N = ^{1 + T3), which is associated to the 
occupation number observable, 

(Nit)) = {^it)\N\ipit)) = {ipiO)uHt,0)\mUit,OMO)), 

(120) 



one recovers eq. (|117p . 

The Grassmann algebra (at a given t) contains only two 
basis elements, {gr} = {gi,go} = {liV"} and the state is 
represented by 

git) = qi{t) + qo{t)^, 

qi{t) = cos(wt)^pi(0) - sin(w<)s^po(0), 

qo{t) = sin(wt)v^pi(0) + cos{ujt)sy^po{0), 

piit) = q,{t)^ , po{t)=qo{t)^. (121) 

The Grassmann evolution operator 

U(t) — cos(ujt) + siniujt) ( — 

= cos(wt) + sin(wt)(a - a^) (122) 
describes the time evolution 

g{t)^U{t)g{S}). (123) 
The conjugate basis elements read {^t} = {i^, 1}, with 

m = qo{t) + qi{t)lp. (124) 



Using (^-^)' 



= — 1 we can write 



u{t) = exp <; w ( V - ^ u 



(125) 



This can be cast into a Hamiltonian form similar to quan- 
tum mechanics, 

idtU = nU , exp(-i-Ht), 

n = iuj(^ij- ^iuj{a-a^), (126) 



with 



0,0 
1,0 



, a' 



0,1 
0,0 



(127) 



We recover the evolution ([95)) . (^7]). demonstrating 

that the functional integral defined by eq. ([7^ indeed de- 
scribes the time evolution of classical probabilities. 



In the Grassmann formulation the occupation number is 
represented by the Grassmann operator 



d_ 

dtp 



(128) 



It translates to the quantum operator N ~ (l + r3)/2, with 



{N{t)) 



d4,~g{t)-^Mt) 



(129) 



At this stage we do not yet have a full description of two- 
state quantum mechanics for the non-commuting spin ob- 
servables in the different directions, even though a non- 
commuting operator appears in the form of % in eq. (|126l) , 
cf. ref. [2] . The implementation of non-commuting observ- 
ables within a classical statistical ensemble, including spin, 
is discussed in detail in ref. ISL 



VIII. GRASSMANN FUNCTIONAL INTEGRAL 
FROM CLASSICAL PROBABILITIES 

The inverse road from classical probabilities {pr} to a 
Grassmann functional integral is also open. In this section 
we start with a classical statistical ensemble for two states 
whose probability distribution follows the time evolution 
(I116p . Using the associated wave function qr{t) and the 
associated Grassmann element g{t) given by eq. (I121[) . we 
realize that this ensemble can be interpreted as two-state 
quantum mechanics. We employ the conjugate Grassmann 
variables in order to derive the associated functional in- 
tegral. This reverses the steps of sect. IVII and demon- 
strates explicitly how a Grassmann functional integral ob- 
tains from a given time evolution of classical probabilities. 
The construction of the functional integral follows the stan- 
dard concept of decomposing the evolution operator into a 
product of infinitesimal factors ^7|. We display it here in 
our formulation in order to show how the functional inte- 
gral arises naturally if we want to express U (t) in eq. (|125p 
in terms of conjugate Grassmann variables. In particular, 
we will not need the often employed concept of coherent 
states. The formalism of this section is easily generalized 
to arbitrary real or complex wave functions for states char- 
acterized by occupation numbers with unitary time 
evolution according to eq. (|123|) . 

In terms of the conjugate Grassmann variables ip the 
basis elements conjugate to gr = {1, ip} a-re given by {gr} = 
{1,V'}, and the Grassmann element conjugate to g{t) = 
qiit) + qo{t)fp reads 



g{t)^qi{t)+qo{t)^p. 



(130) 



We want to implement the time evolution operator (jl25l) in 
the formulation with conjugate Grassmann variables ipit). 
In order to perform the necessary normal ordering of the 
factors ip and d/dip in the exponential (|125p we proceed by 
the standard methods of a formulation of a path integral. 
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This will bring us back to the functional integral ([83]) and 
generalizations thereof. The central idea is a decomposition 
of U{t,tQ) into a product of operators with infinitesimal 
time steps U{tn+i,tn), tn+i — i,i = e — ^ 0, such that the 
appropriate ordering for U{tn+i,tn) becomes trivial up to 
negligible corrections in higher orders in e. 

1. Factorization of evolution operator 

We first proceed to a split of the evolution operator into 
factors for subintervals in time. These intervals can be 
taken infinitesimally small. For an arbitrary operator 
■Flip] = a + + c-^ + d^ip and an arbitrary element 
5'(V') — 9 ^ hip one has the identity 



di/. e^'l' T[i>]g{ip) (131) 



with T[ip'] = a + btp' + c-^ + d-^ip' . This can be used 
for a decomposition of the product of two operators 



diPe^^T,[^,]T2m{iP) (132) 
dij'e'^^'Tiiij'] f dVrf^'e^'^^~^'^^2M5(V')- 



(The product is here the standard product of Grassmann 
operators and should not be confounded with the product 
o used in eqs. (11141) . (jllSp .) In particular, we can employ 
(for arbitrary i8['0]) 

BU{t,to) = , Ti^BU{t,ti), 

T2 = U(tiM). (133) 

and identify ip' — ?/;(io),'0' = V'(^i):V' = ■0(^o), such that 
dip S'^{m{tM))W\9{^) = I diP{h)e^^^''^ 

{BU{t,h))mh)] J dV(to)dV'(io)e-'^^*")('^(*^^-'^(*"^) 
U{tutomto)]9ii^(to)). (134) 

For ti — to = e the integral / dip (to) dip (to) involves in the 
integrand only an infinitesimal evolution operator that can 
be evaluated in lowest non-trivial order in e. 

This procedure can be repeated for BU{t,ti). We end 
with with a sequence of simple Grassmann integrals over 
infinitesimal evolution operators 



dip e^^{BUltM)[4']9{4>) = J dm e'^^^'^Bim] 

N-l „ 

n { / diPHMitn) e-^(*")(^(*'^+^)-^(*")) 

n )]}9{m))- (135) 



Here = t and factors with larger t„ are to the left of 
factors with smaller <„ . We recognize in the first factor the 
"shift operator" 

S{t^+M = |dV'(i„)e-'''(*")(^'*"-^^-^(*")) 

= d{m)-Htn+l)) (136) 

with multiplication law 

S{tn+l,tn)og{lP{tn))= (137) 
dlp{tn)S{tn+l,tn)9{lp{tn)) = g{lp{tn+l)). 



Its role is to shift the argument of and we recognize in 
eq. (I135P the product o with factors 5(t„+i, t„)W(t„+i, i„), 
which first apply the Grassmann evolution operator U and 
subsequently shift the argument. For B = 1 this yields the 
factorization of U, 



diPito)e'^^^''>^Uit,to)[iPito)]g{iPito))^ J diP{t)e^^^'^ 
) oU{t„+i,tn)[m)])g{iP{to)). 



n=0 



(138) 



2. Hamilton operator 

In the limit e — tn+i — in it is sufficient to consider 
the lowest order in e, 

U{tr,+i,tn)mtn)] = 1 - in[d/di/j{tn),ip{tn)]e. (139) 

Since H is multiplied by a ^-function from S^(t„+i,t„) we 
can replace the argument ip(tn) by ip{tn+i). Furthermore, 
if H is normal ordered, as for eq. (|126p . we can replace 
d/dijj{t„) by ip{tn). In total, we perform the replacement 

U{t^+i,t^)[iP{t^)] ^ 1 - ii?[^(t„), V(i„+i)]e 

= exp| -ii/[V3(i„),^(i„+i)]e}. (140) 

For the particular example of two quantum mechanics this 
yields 

U{tn+i,tn) exp[ea;(V'(i„+i) - ip{tn)]- (141) 
We can now use eq. (|138p for a computation of g{t), 

J diP{to)e^^^'MtM)[m)\9{to) 

= j diP{t)e^^^'^9{t). (142) 

With 



U{tn+l,tn)g{lp{tn)) = / dlp{tn)u{tn+l,tn)g{lp{tn)) 



(143) 
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and 

one obtains 
9{t) 



exp 



d^jj{tn) exp I - ieH(il}{tn), ip{tn+i) I 
{^(i„)(V'(i«)-^(i„+i))} (144) 



Z^(t,io)[^(io)]5(io) 



(145) 



Here ^(t) depends on ^p{t) and g(io) depends on ip{to), 
while multiplication with U{tn+i,tn) involves a shift from 
il){tn) to ijj{tn+i) where the product o is employed. We 
recognize eq. (|115l) . The derivation of eq. (|142p has not 
used any specific properties of lA except the existence of 
an expansion (|139l) in lowest order in e. The results of 
this section can be used for arbitrary U or %, with trivial 
generalization to g constructed from an arbitrary number 
of Grassmann variables. 

3. Functional integral for observables 

We want to express the expectation value of some observ- 
able {A{t)) in terms of the initial wave function f;(-0(to)). 
According to eqs. p^ . ([TB| we write 



{A{t)) = / dmmAglt) 

d^{to)g{to)U{t, to)AU{t, to)g{to), (146) 
and express with eq. pO)) 

-gito) = / d#(^)e'^'^(*°). (147) 



We therefore need an expression for l({t,to)A{t)l4{t,to)- 
For this purpose we consider in eq. (jl33l) 



Bm)] = {u{to,t)A)[m] 



(148) 



We can continue the factorization procedure and end with 
the main result of this section 

{A{t)) = I d^^drgi^) J'^{U{tM)AU{tM)Mg{^) 



oV'(*2n)V'(*2n) 



W [ / dil){tk)di^{tk)\g{'4'{t2N)) 

N-l 

m ), V'(^m+l)](im — ^m+l)} 



k=0 
2JV-1 

X 

=N 



X exp { - 'ilj{tm){i^{t 



m+1) 



X Yl [exp { - ii?['0(i„), ■!/'(t„+i)](i„+i - t„)} 

Xexp{-i^{tn){^{tn+l)-^{tn))}]g{i^{to)), (149) 



where t = In and t/j = iip{t2N)- The order of the factors is 
with increasing n and m from right to left. Counting the 
t„i backwards, t,n+i — tm = — e we associate t2N — to- For 
a given t„ = io + we end then with two distinct Grass- 
mann variables 'ipitn) and (p{tn) — i^{t2N-n) (except for 
'0(^-/v)-) This amounts to the Schwinger-Keldish formalism 
in quantum field theory f9|. 

The structure of eq. (|149p can also be understood from 
the equivalent expression 



{Ait)) - / dmdme*^'''^^'^mA[m,m]g{t), 

(150) 



where g{t) is given by eq. (jl45p and g{t) obeys 



m=g(.to) n (°^^(^n,t„+l)). 



(151) 



In this expression g{t) depends on ifi{t) — 4'it) and 5(^0) on 
If (to), while multiplication from the right with V{tn,tn+i) 
involves a shift from arguments (p{tn) to ip{tn+i)- The 
order of factors in eq. (|151l) is now with smaller t„ to the 
left. The product involves an integration over ip{tn), 



f{tp{tn)) O V{tn,tn+l) = J dip{tn) f {(p{tn)) v{tn, t^+l) ■ 

(152) 

For a hermitean Hamiltonian one has 

v{tn, tn+l) = v{(p{tn), (p{tn+l)) 

d(p{tn) exp ^^(p{tn){(pitn) - <^(i„+l))| 

exp|ieiJ((^(t„+i),^(t„))|. (153) 

Here / obtains from / by changing the sign of all odd 
products of Grassman variables. 

The fact that U{t, to) is an inverse oiU{t, to) is refiected 
by the relation (for (p{t) = ip{t)) 



dil;{t)d-^{t)e^^'^^'^'^ j d^{t-e)f{ip{t-e))v{t~e,t) 
X I d4^{t~e)u{t,t-e)g{^/j{t~e)) 

diPit -e)J difiit - e)e^(*-^)'/'(*-0 
xf{0{t-e))g{^{t~e)). (154) 

For arbitrary t' < t and arbitrary g{t') = g(^jj{t')) and 
fit') = f{ip{t')) one has 

diA(i)#(t)e^(*)^(*'/(i') 



t-e 



W V{tn',tn- + l) n U{tn+l:tn)og{t') 
t„,=t' t„=t' 



dil;{t')d(p{t')e'^ 



^f{t')git'). 



(155) 
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A convenient way to understand these relations relies on 
the familiar concept of transition amplitudes which can be 
defined in the Grassmann formalism as 

(/|5> = / dHi^e^^'fWgW - J (156) 

with {g\g) = 1. In this language the relation (|154p reads 

ifitMt)) = {f{t^e)oV{t-e,t)\U{t,t-e)og{t-e)) 
= iht-eMt-e)), (157) 

while eq. (|155p can be expressed as 

ifitMt)) ^{fit'Mt')). (158) 

For a Grassmann operator J-]^id/dijj,ip) we generalize eq. 

m 

ih^Ng) = iI\^N\g) 

= J Di;fi^)FNii^,^)gi^), (159) 

where it is understood that d/dtp - and ■0- factors in 
J^Ar(9/9V', ip) and FnH^, ip) are to the left of -i/'-factors. 

4. Functional integral with future time 

The form of eq. (|149p is suggesting a different option, 
namely to count t^ to the future, t^+i — tm = ^, asso- 
ciating t2N — 2t — tQ. We will see below that the cor- 
responding change of the sign of e in the terms involving 
H[ipitm),ipitm+i)] can be absorbed by a simple change 
of g('0(i2Af))- Instead of being conjugate to gito) it is now 
conjugate to g(t2Af) = gi'^t — to). In consequence, eq. (|149p 
yields the functional integral expression 

iAit)) = J VijViPit)git + s)Am),m] 

xf{exp{-^i(t')}}5(t-s), (160) 
t' 

with git + s) depending on ipit + s) and g{t — s) involving 
ipit — s), whereby 

Lit') ^ i^it'){i:it'+e)-^Pit'))+teH[4^it'),i:it'+e)]. (161) 

In cq. (|160p wc integrate over all ipit'), "tpit') and t' covers 
2N + 1 values in the interval t — s<t'<t + s, s = Ne, 
with Lit + s) — —tpit + s)'!/'(i + s). We have assumed that 
A is a bosonic observable, containing only terms with an 
even number of ipjip- The time ordering T is only needed 
if H contains terms with odd powers of ip, ip, such that 
H{ta) and Hitb) do not commute. The action of T is to 
put factors exp— iei?(t') with larger t'-arguments to the 
left of factors with smaller t' arguments. 

We recall that g{t + s) is conjugate to git + s) and g{t + s) 
is related to git — s) by the evolution law. Therefore only 
one independent element of the Grassmann algebra spec- 
ifics the state of the system. The expression (|160p equals 



the functional integral ([83l) (for Z = 1) if we identify 
gin = git — s) and cjf = git + s). This direct construc- 
tion of the relation between gf and gin coincides with the 
requirements for a consistent evolution of git) and git) de- 
fined by eq. (|88)) and discussed in sect. |Vll For our specific 
example we have derived the Grassmann representation for 
the simple time evolution for classical probabilities (|116l) . 
The formulation in terms of "0 and '0 leads directly to the 
functional integral expression (1551) in sect. I VII 

What remains to be shown is that the functional integral 
(I149P with^m+i > tm and5(0(i2jv) = git+s) indeed yields 
the expression (jl46l) or (I150p . For this procedure we write 
eq. (jl50p in a form that can directly be used in eq. (jl55p 

it' — tN,t — t2N), 

iA{t)) = J dVXtjv)#(iiv)e^(*")'^(*")g(tjv)(-4.9)(ijv) 
= J dijit2N)d4'it2N)e^^''-^^^''-^gitN) 

i2N-l t2N-l 

X n yitn',tn' + l) n U{tn+l,tn)iAg)itN) (162) 

With 

t2N-l 

gitN) n Vitn',tn,+i)^git2N), (163) 

and replacing again A[^/^'^p{tN),^lJit]y)^ 
A[0(i7v), '0(iAr)] , we find indeed eq. (|149p with t^+i > tm 
and gi^pit2N) = git2N) = g{t + s). 

IX. TWO COMPONENT SPINOR AND 
COMPLEX STRUCTURE 

Our first example has involved a real wave function and 
a real Grassmann algebra. In this section we consider a 
simple system that admits a complex structure. As a con- 
sequence, we can equivalently use a complex wave function 
and a complex Grassmann algebra for a description of this 
system. The opposite mapping from the complex descrip- 
tion to the real one is rather trivial, since complex numbers 
can always be expressed in terms of a pair of real numbers. 

1. Four state quantum mechanics 

Our example describes a four state system, r — 1 . . . Ms , 
Afs =4, B = 2. The time evolution of the classical proba- 
bilities is given by 

Piit) = Pi,o, 

P2(i) = coa^ {ujt)p2o + siu^ i(^t)p3o 
-2 cos(wt) sin(a;t)^p20P30, 

Psit) = cos^iu}t)p3o + sin^{ujt)p2o 
+2 cos(wt) sin(wt)^p20P30, 

Piit) = P40, (164) 
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where Pr,o are the initial probabihties at t = 0, normalized 
by P10+P20+P30+P40 — 1- This evolution can be described 
by a differential equation for the real wave function with 
components {^rjiPr — <1^, 

dtqi = dtqi = , dtq2 = -^93 , dtqs = i^q2- (165) 



We observe the close analogy with sect. IVIl with (72 and ^3 
playing the role of qi and go in eq. (|119p . 

However, we have now a four-state system with other 
possible observables. It admits a Grassmann representa- 
tion different from sect. IVIl A suitable basis for a real 
Grassmann algebra consists of {gr} = (1, V'l? '02, '0iV'2)- In 
this basis the Hamiltonian for our example can be written 
as 



d d Q 

'H = -iuJ^T—1p2 + ZW— — 01 = w TT— (t2)q/3'0/3- (166) 



Expressed in terms of the conjugate Grassmann variables, 
in a notation with two component spinors 0, it reads 



H = UJIp T20. 



(167) 



This can directly be used for the construction of the 
corresponding functional integral, with tp taken at i -I- e 
and ?/; at t as in eq. ([M|) . 



(168) 



We note that H is now a bosonic quantity such that no 
time ordering T is needed for the functional integral (1771) . 
with action 



S = Y. {^(tfi^i* + - 0(0) + zeif [0(i),0(i + e)] }. 



(169) 



As it should be, the Hamiltonian (|167p is hermitean 



(//[0,0])^ =i7[0,V]. 



(170) 



In the continuum limit we can write the action as a sum 
of a "dynamical term" Sdyn and a "Hamiltonian term" Sh, 



S = Sdyn + Sn , Sdyn^ j dtiPitfOtm, (171) 
Sh = i ( dtH[t) , H{t) = H{'ijj{t),'il;{t)) ^ H\t). 



This implies antihermiticity of S — —S\ and therefore 
hermiticity of the Minkowski action (|76l) . So far the action 
S of our example is real and we therefore consider a real 
Grassmann algebra. 

2. Complex structure 



The time evolution (|165p can also be described in a complex 
basis 



(Up to normalization, this corresponds to a spin in a con- 
stant magnetic field.) The definition of the complex two- 
component wave function q in eq. ()172|) defines a complex 
structure for the four real components qr- It is given by an 
involution 0, whereby 93 and 54 are odd with respect to the 
involution (they change sign), whereas qi and q2 are even 
(invariant). Not all possible observables for the real wave 
function {^t} are compatible with the complex structure. 
Compatibility means that observables can be expressed as 
complex operators A acting on without involving the 
complex conjugate q* . In other words, the expectation 
value of compatible observables obeys 



{A) = q^Aq 



(173) 



An example for an observable which is compatible with the 
complex structure (|172l) is given in the real basis by {Ar} — 
{0, 1,1,0}. In the quantum formalism this corresponds to 
the real diagonal operator 



A = diag{0, 1, 1,0) = (A^i - iVz)^ 



(174) 



It equals one if iVi and N2 are different, and zero if they 
are the same. In the complex basis it is represented as 
A — (1 — T3)/2. Since H — —loA the quantum formalism 
implies that {A) is time independent. This can, of course, 
also be inferred directly from the time evolution of the 
probability distribution (I164p and the values A^ for the 
four states of the classical ensemble. 

The complex structure for the wave function is reflected 
in a possible complex structure for the Grassmann algebra. 
The two-component wave function (gi, 92) = (^i +«94, Q2 + 
iqz) can now be associated to an element of a complex 
Grassmann algebra 



g = qi+ 920- 



(175) 



This algebra uses only one independent Grassmann vari- 
able instead of the two variables 0i , 02 for the formula- 
tion in terms of a real Grassmann algebra. In the complex 
formulation the Schrodinger equation takes the form 



idtg = Ug = q2ip, 
d 



(176) 

H - (^^0-lj =c^(AA-l). (177) 

In terms of the conjugate Grassmann variable this yields 

H = w(i/;^-l), H^=H, (178) 
and we can formulate the functional integral using 



Sh = ii^ / dt(0(t)0(O - 1). 



(179) 



All steps of sects. IIVI and rVTl remain valid for this complex 
Grassmann algebra. 

The equivalence of a complex wave function and Grass- 
mann algebra with a real wave function and Grassmann 
algebra is easily generalized 0. The real formulation has 
one Grassmann variable more than the complex one. Out 
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of the B bits in the real formulation one bit can be used 
for differentiating between the real and imaginary parts of 
a corresponding complex wave function. The remaining 
B — 1 bits can be used for the construction of a basis of 
a 2^~^-component complex wave function. The latter is 
equivalent to a 2^-component real wave function. In con- 
sequence, real wave functions {qt} can cover the most gen- 
eral case of quantum wave functions. The complex wave 
functions and complex Grassmann algebra which are char- 
acteristic for many quantum systems reflect the presence 
of an "additional" complex structure. 

We finally observe that an analytic continuation to "eu- 
clidean time" t = it leaves Sdyn invariant, while Sh picks 
up a factor —i from dt (assuming Sh contains no time 
derivatives) 



Se 



dT{i;dri^ + H{4>,^j)). 



(180) 



The corresponding euclidean functional integral realizes 
Osterwalder-Schrader positivity [l^l ■ 



X. RELATING PAST AND FUTURE 

The formulation of a functional integral with one 
"boundary condition" set in the past, and the other 
one g / set in the future needs the relation between g f and 
gin- If the measured observables and correlations involve 
only observables within a certain time interval, we can use 
gin {tin ),gf{tf) with tin and tf at the boundaries of this 
interval. If we start from a functional integral which in- 
volves arbitrary times {tin — ^ —oo]tf — >■ oo), the wave 
functions gin{tin) and gf{tf) are the only pieces of infor- 
mation needed from the probability distribution outside 
the interval which is relevant for the observables [l^l- Still, 
for the most general case the relation between gf and gin 
can be complicated, which would make a practical use of 
this formulation difficult. 

1. Time translations 

For many relevant systems, however, the dynamics is char- 
acterized by a time independent Hamiltonian H or evolu- 
tion operator K.. In the functional integral formalism this 
corresponds to an action that is invariant under time trans- 
lations. In this case the factor g/ is easily computable for 
a large class of states specified by the initial wave function 
gin- If we assume the presence of a complex structure it is 
sufficient that gin is an eigenvalue of the Hamiltonian 



"Hgin — Egin- 
This implies the simple time evolution 



Cr{t) 



-iE{t-tir^) 



(181) 



(182) 



Thus Cr{tf),g{tf) and g{tf) are easily obtained from gin 
and the formulation of the functional integral (jl60l) can be 
employed without major obstacles in practice. A particu- 
larly simple example is the computation of vacuum proper- 
ties, for which g is time independent and therefore ^/ = gin- 



One also can choose tf — tin = jE — t such that 5/ = 
zte'^'^^gin,gf = ±e~*'^^gi„, with the minus sign realized for 
m odd. In this case one can put time on a torus, with 
periodic or antiperiodic boundary conditions on g{t) for m 
even or odd. Indeed, we can write the functional integral 
([T601) as 

{A{t))= fv^bV^pA[i;{t),^{t)]e-^'p{t,n), (183) 



where we assume for simplicity that A and S involve only 
even powers of Grassmann variables. We identify tf = tin — 
e such that the functional integral is periodic in time with 
period mn/E. We note that p{tin) only involves Cr{tin) but 
depends on the Grassmann variables ip{tf) = '4'{Un — () and 
'4'{tin)- Furthermore, we want S to be translation invariant 
around the torus such that we have to supplement the miss- 
ing terms exp | - (<i„ - t)ip{tin )-ieH[tp {tin - e) , "0 {Un)] | ■ 
We compensate this by the definition of p{tin), 

p{tin) = (-l)™^C*(tm)cT(tm)exp{Vj(ii„ - e)'(/'(im)} 

xexp^^ie(^H[ip{tra - e),1p{t^n)] -e)} 
y-9p{i^{Un-e))gr{i^{tui)), (184) 
with cjp conjugate to gp, such that 

dip{tin)dip{tin - e) exp {%l){tin - e)^{tin)] 

xgp{^{Un - e))gr('0(ijn)) = Srp, (185) 

and / dlP{t,n)di:{Un - €)p{t,n) = ("l)". 

The expression (|183p computes {A{t)) as a correlation 
function for the "observables" A and p{tin) in a periodic 
functional integral. All information about the particular 
state is encoded in p{tin)- We can easily generalize eq. 
(|183l) to observables A involving several time arguments ti 
inside the interval tin < ti < tin + rmr/E- This type of ob- 
servables for energy eigenstates (|181l) can be obtained with 
periodic time, and a "long time behavior" with At > rmr/E 
never matters. In analogy to the Matsubara formalism one 
may switch to frequency space with discrete frequencies 
ujn = (2n + l)E/m,n G TL- We recall, however, that one 
works now with "real time" without analytic continuation. 

The presence of a complex structure is quite generic in 
our setting. Besides structures of the type discussed in the 
preceeding section, we can use charge conjugation as the 
involution underlying the complex structure. For our Ising 
type classical statistical systems charge conjugation maps 
empty bits to occupied bits and vice versa, r = [ric] 
f = [na] = [(1 — Jt-q)] • On the level of the real Grassmann 
algebra we can define the charge conjugation as 



1 , = 1. 



(186) 



This extends to arbitrary basis elements of the Grassmann 
algebra. 



(187) 
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for an appropriate choice of signs ar = il- If the evolution 
/C commutes with C we can use C for the definition of a 
complex structure - an example will be given in the next 
section. 

2. Time reversal 

The elements and gin{tin) may also be related by 

time reflection. In general, the dynamical term in the 
action is invariant under a time reflection around tg = 

(t/+i™)/2, 

/3 

l/'a(t.+i)^^ra/3^/j(is-0, (188) 

/3 

where T and T are related by T"^ = ^ accompanied by 
a total reordering of all Grassmann variables. The time 
reflection is an involution, = 1. If for an appropriate 
choice of T and T also the Hamiltonian part Sh remains 
invariant, the action is time reflection invariant T{S) — T. 
For the action (|74|) this is realized by T = T = —1, while 
for the action (|169p we use in eq. (|188|) the matrices 

T^f^Ti. (189) 

(Om deflnitions are such that reflection symmetry holds in 
the formulation with discrete time steps t„ without further 
modifications.) The functional integral is T-invariant if 

T{gf{tf)) = ±gUUn) , r(5m(tm)) = ±9f{tf)- (190) 

We emphasize that time reflection symmetry is a property 
of the whole family of probability distributions {pr(^)} lo'^ 
all t, rather than being a symmetry of a state at a given t. 

For a time reflection invariant action the condition (|190p 
is obeyed if 

r(5(t,)) = ±g%). (191) 

This follows from the reflection symmetry of the time evo- 
lution. The easiest way to see that eq. (|19ip follows from 
eq. (|190p relies on the definition (|90p . realizing that S'< is 
mapped by time reflection to 5>,T(e~'^<) = e~^> . The 
inverse follows from the invertibility of the time evolution. 
For our flrst example and is = a time reflection sym- 
metric setting is realized for g(Q) = 1 or 5(0) = tp, cor- 
responding to pi(0) = 1 or po(0) — 1. One verifles that 
the corresponding time evolution of the probability distri- 
bution (|116p is indeed reflection symmetric. For the sec- 
ond example in eq. IIXI we have reflection symmetry for 
ff(0) = ^(?/'i ± i/'2), corresponding to p2{t) = Pz{t) in eq. 
(|164l) . If eq. (|19ip is obeyed for some suitable ts we may 
choose tf and tin obeying tg — {tf + ii„)/2. Then eq. 
()190p allows us to compute gf{tf) in terms oi giniUn) and 
the functional integral is fully specified by the coefficients 
Cr{tin)- We observe in this context that the choice of the 
matrices T and T in eq. (|188p . which leave S reflection 
invariant, may not be unique. Any choice which realizes 
eq. pMt wfll do. 



XI. QUANTUM FIELD THEORY FOR 
TWO-DIMENSIONAL FERMIONS 

In this section we generalize our previous examples to a 
quantum held theory for free fermions in two dimensions 
(one space and one time dimension). We will construct a 
map between the functional integral defining a quantum 
field theory for fermions and an associated probability dis- 
tribution of a classical statistical ensemble for Ising type 
variables. Further generalizations to fermionic quantum 
field theories with interactions and in arbitrary dimensions 
are straightforward. We use a real Grassmann algebra and 
the index a covers a space index Xi on a one-dimensional 
lattice with x^+i — Xi — e, as well as an internal index 
7 = ± . We may use a torus with L points Xi and circumfer- 
ence I = Le, such that B = 2L. In the continuum limit we 
deal with two Grassmann valued fields ilj+{t,x),ijj-{t,x). 
The functional integral involves only real elements, even 
though we use i occasionally for the purpose of analogy 
with quantum mechanics. 

1. Functional integral and Schrodinger equation 

In a discrete setting the action in the functional integral 
formulation is given by 

^ = (^+(^ + e,a; - e) - tp+{t,x)) 

t,X 

+ tp- {t, x-e) {ip^ [t + e,x)- tp- {t, x-e))} 

(192) 

or equivalently 

S = x) (^+ {t + e,x)- ij+ {t, x)) 

t.x 

-I- i/)- (t, x-e) (■(/'_ {t + e,x - f) - ip- {t, x-e)) 

- 4'+{t, x) {tp+{t + e,x)-ip+{t + e,x-e)) 

+ V'_(t, x-e) {ip^{t + e,x)-'tp^{t + e,x- e)) }. 

(193) 

The continuum limit involves now a rescaling of the Grass- 
mann variables ip±{t,x) — >■ y/eip±{t,x) such that 

S = {ip+dt^+ + ip-dtip^ - il'+dxilJ+ + ^-dxip~ ] 

Jt,x 

= f ^P^dt^p + 1 ( H, (194) 

with ip — (-04- ,?/>„) , f/;^ = (t/i 4-, -!/;_). The Hamiltonian part 
reads 

H = i j {'^+d^i>+ - i>-d^^-] = i f VVsa.V, (195) 

and obeys W ^ H, H* = = -H, such that = -S 
and S\j — Sm ■ The functional integral is defined as in the 
preceding sections. 
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The Schrodinger equation for the Grassmann wave func- 
tion g obtains as 

dtg = -iUg, 



where we replace V'i in eq. (|195p by d/dijj±. The Hamilton 
operator T-L commutes with the two "particle numbers" 



(197) 



Thus A/+ and AC are conserved. We can consider simul- 
taneous eigenvectors of A/+ and A/"- 

(198) 

and expand a general wave function as 



9 — ^ ^ -^n-i- ,n_ ^71+ ,n_ • 
,n_ 



(199) 



2. Multiparticle states 

The one particle state with = 1 , n_ = can be written 
as 



51,0 = / <li,o{x) 



d 



dtp+{x) 



with vacuum state 



|0) - / qi.^o{x)aUx)\0) (200) 



(201) 



= n^+^^^'^-^^^- 



Here we also employ the notation of sect. IIIII by associ- 
ating d/dip with the creation operator . The real one 
particle wave function is given by 91,0(^1 x) and normalized 
according to 



9i,o(^) = 1- 



(202) 



On the torus we use antiperiodic boundary conditions for 
the wave functions as qi^o{x), such that the associated 
probabilities as pi^oi^) are periodic. 

The wave function gi^o obeys a one-particle Schrodinger 
equation which is derived by inserting g — gi Q in eq. (|196p . 
This becomes a very simple differential equation 



dtqifi{t,x) = dxqi.o{t,x). 



(203) 



The general solution (71,0(2; -f t) depends only on x -I- 1 and 
therefore describes a "left-moving" particle. The arbitrary 
shape of the wave function is conserved and specified at 
some initial time by qi^Q{tQ,x). Similarly, 50,1 describes 



"right-moving" particles for which the real wave function 
'Zo,i(a; — t), which solves the Schrodinger equation 



(204) 



depends only on x — t. 

Two left-moving particles involve the two-particle wave 
function q2fi{x,y), 

\nJx..y dip+{x) d-ip+iy) 

with normalization 



(206) 



As appropriate for two identical fermions it is antisymmet- 
ric under particle exchange 



(207) 



In contrast, one left-moving particle at x and one right- 
moving particle at y differs from the state where the left- 
mover is at y and the right-mover at x. Thus for 



51,1 = 



d d 

gi.i(a:,2/)— — -— — -|0) (208) 

di^+{x) dij^{y) 



the wave functions qi^i{x,y) and qi^i{y,x) are not related 
by symmetry. We encounter the standard situation for 
different particle species. 

The Schrodinger equation for two left-movers becomes 

dtq2fi{t, X, y) ^ {dx + dy)q2fi{t, x, y). (209) 

A special solution is given by a factorized wave function 

q2,o{t, X, y) a{x + t)h{y + t)- a{y + t)b{x + t). (210) 

More general solutions obtain by superposition of terms 
of the type (|210p . The generalization to states with n+ 
left-movers and n_ right-movers is straightforward. Our 
normalizations are chosen such that the general Grassmann 
element ()199p obeys 



(211) 



It is obvious that qn^,o is the wave function for n+ identi- 
cal fermions, with total antisymmetry with respect to the 
exchange of any two particles. 

3. Classical wave function and probability 
distribution 

We next turn to the classical statistical ensemble which 
is described by the action (|192p . with probability distri- 
bution {pr} and associated classical wave function {qr}- 
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The general wave function qrit) evolves according to the 
Schrodinger equation obtained from eq. (|196p by inserting 



lit) =^<}r{t)gr 



(212) 



The states r of the classical statistical ensemble are 
given by the sequence of occupation numbers, r — 
(x) , n_ {x)] . The particle numbers , n_ in a given 
state T obey 



For any given time t we decompose 



(213) 



(214) 



where gn+,n^ is normalized such that Pn+,n_ ~ ^n+,n_ 
equals the probability to find n+ left-moving particles and 
n_ right- moving particles. Since the particle numbers Nj^^ 
iV_ are conserved, the solutions of the Schrodinger equa- 
tion have time- independent amplitudes In turn, 
the general solution for the fixed particle wave function 
qn+,n-, defined by 

5n+,«_(i) = (n+!n_!)-^ 

X / <7n+,n_(a;i . . .a;„_|_,2/i . . .2/„_,<) 

Jxi...x„^,yi...y„_ 

(215) 



with 



d 



d 



d 



d 



(216) 

|0) 

and 



have been discussed above. The coefficients An 
the wave functions qn+,n^(t) fix completely the evolution 
of the probability distribution Prit). 

In the opposite direction we may start with a given solu- 
tion pr (t) of the evolution equation. From this probability 
distribution we can compute directly the conserved prob- 
abilities Pn+,n- by summing the probabilities of all states 
with a given number n+ of occupied bits of type -I- and 
similar for the type — . Without loss of generality we can 
choose a convention where A„^^„_ — ^Pn+,n^ ^ 0. From 
the relative probabilities to find the n+ occupied bits at po- 
sitions Xi . . . Xn+ and the n_ occupied bits at yi . . . ?/„_ we 
can determine the wave function qn^,n- [xi ■ ■ ■ yn- ) up to a 
sign. The overall signs of the wave functions qn+,n^{ta) at 
some initial time to are free. The expectation values and 
correlations of observables do not depend on these signs 
and we may choose them according to some arbitrary con- 
vention. (This is similar to a choice of gauge). For all times 
t the sign of qn+,n- (t) is then related to (/„_,.,„_ (to) by the 
solution of the Schrodinger equation. Furthermore, at to 
only the overall sign of 9n+,n_(^o) is arbitrary. The rel- 
ative signs for different arguments (xi, . . . , yi, . . . y„_ ) 



can be fixed by requirements of continuity and differentia- 
bility (3|. Up to some irrelevant "choice of gauge" the wave 
function {qrit)} for a system of free particles in two dimen- 
sions is determined by the probability distribution {prit)} 
of a classical statistical ensemble. For our example, the real 
quantum wave function ipgit) = J2n+,n- ^n+,n_gn+,n_ (0 
for a multi-fermion system is directly related to the classi- 
cal wave function {qr}- If the time evolution of the prob- 
ability distribution {pr} is such that the associated Grass- 
mann element (I214p obeys eq. (|196p . we can express all 
relevant expectation values by a functional integral with 
action (fTM]) . 

4. Lorentz symmetry 

Our system of free fermions is invariant with respect to two- 
dimensional Lorentz transformations. This is most easily 
seen by using 



In terms of the real Dirac matrices 

7" = «'^2, 71 = n, 

which obey 

{7^7''} = 2ry^^ 
with r]^'^ = diag(—\, 1) we can write 

S = - t V^7^9^^. 



(217) 
(218) 
(219) 

(220) 



The infinitesimal Lorentz-transformations act on the real 
two component spinors ip and '0 as 



with 



S" = -i[/,7^]=-^r3. 



01 



(221) 



(222) 



They commute with 7 obeying 

7 = 7V = r3, 7'-!, {7^7}- 0. (223) 

The eigenstates of 7 correspond therefore to irreducible 
representations of the Lorentz group 



(0,0;+) = ^^(1-7), (-^-,0) - ^7^(1 + 7). 

(224) 

The Weyl-spinors i/'i transform by multiplicative rescal- 
ings 



H+ = 1-0+, Si:^ = 



(225) 
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In addition, the coordinates 
the usual way. 



(t, x) are transformed in 



5. Discrete symmetries 

Besides the invariance under the continuous Lorentz- 
transformations the action (|194p . (|220l) is also invariant 
under the discrete symmetries of parity reflection (P), time 
reversal (T) and charge conjugation (C). Parity reflection 
acts by reversing the sign of the space coordinate, x — ^ —x, 
accompanied by a discrete map acting on the Grassmann 
variables 

P {i,{t, x)) = -7 -x) = rMt, -x), 

,- . - „ _ (22b) 

P{ij{t,x)) =-tP{t,-x)jj'' = -ij{t,-x)T,. 

It exchanges two Weyl spinors V'+ , ip- , 

P : tP+^iP^, (227) 
Time reversal changes the sign of t and maps 

T: V_f^7/>+. (228) 

In addition, all Grassmann variables are totally reordered 
(transposition). 

Gharge conjugation exchanges the Grassmann variables 
with their conjugate ones (no action on the coordinates) 



(229) 



— ^,0 



(This corresponds to a charge conjugation matrix Ci 
in the general setting of ref. [11[. The use of the same 
name of charge conjugation should not hide the fact that 
the transformation (|229p is conceptually different from the 
charge conjugation discussed in sects. |V] and [X] It acts 
within the extended Grassmann algebra constructed from 
■0 and -0, in contrast to eq. ([57]) .) 

The combined symmetry PTC maps ^-y{t,x) — >■ 
ip'ji—t, — x), 07 (t, x) — '(l:~f{—t, —x), together with a trans- 
position. Furthermore, the action (jl94p is invariant under 
the reflections 



^4 



0_ 



-0+, 
-0- 



(230) 



This allows us to define modified discrete transformations 
as P- = S-P etc. Of course, the action remains invariant 
under these modified transformations as well. Finally, the 
transformation — > 7-0, — tp^, combined with a trans- 
position, leaves S invariant. This may be used for defining 
a version of time reflection without transposition. 

A simple extension of our model could use a complex 
Grassmann algebra with a complex wave function {c,-} in- 
stead of the real {qr}- For free fermions the real and imag- 
inary parts of {cx} evolve independently and the extension 
accounts for a trivial doubling of the degrees of freedom. 
We could now add a pointlikc interaction, S S + Sj, 

Si — iX ip+(t, x)ij)- (i, x)ip- {t + e, x)ijj+{t + e, x) 



iX / '0+V'_0„0'4 



This interaction is invariant under Lorentz-transformations 
as well as the discrete transformations P, T and C. The 
interaction part of the Hamiltonian 



H — X 0+'0_-0_'04 



(232) 



is hermitean only for real A. This is the reason why a 
pointlike interaction is not possible for a real Grassmann 
algebra. In the following we will not discuss this extension. 
We will rather discuss a complex structure which allows to 
represent the real wave function {qr} in terms of a complex 
wave function {cr} with half the number of components. 

6. Weyl spinors 

We observe that the action for free fermions can be written 
as a sum 



s = s+ + s. 



(233) 



where S+ only involves the left-movers "0+ and "0+ , whereas 
S- accounts for the right-movers ip- and ip-. The func- 
tional measure can be written as a product of integrals, 
one involving 04., 0+ and the other 0_, ip^. If the initial 
state factor has a factorized form 

gin = 9 (0(tm)) = 9+ (■0+(tm)) 9- H^-iUn)) (234) 

this also holds for g/ — g {i>{tf)^ and the partition func- 
tion ([77]) factorizes 



z = Z+Z^ 



(235) 



In this case our setting describes independent "worlds" for 
the left-movers and right-movers. The wave function fac- 
torizes 



qr{t) = q^+'^(t)qi-^{t), T ^ {p,a) 



(236) 



with p — [n^{xi)\, a = [n^(xi)]. This extends to the 
probability density 



p.(t)=p(+)(<)p(-)(<), p(±)(i)-(9f^)' 



(237) 



In fact, a model involving only the right-moving Weyl 
spinors 0_, 0_ is fully consistent. It obtains from our 
setting by removing 0+ and 0+, i.e. by putting 0+ = 
■0+ = in the action and omitting 0-|-, ip^ in the func- 
tional measure. The action is then invariant under Lorentz- 
transformations and charge conjugation, but violates P 
and T since these symmetries map the right-movers to the 
left-movers. A model involving only the Weyl spinors 0_, 
ip- still conserves the discrete symmetry 



PT 



0_ O 0_ 



(238) 



(231) 



and is therefore also invariant under PTC. 

In two-dimensions Maj or ana- Weyl spinors are compati- 
ble with Lorentz-transformations |12l |. They obey C0 = 
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and obtain by identifying ■0+ and as well as ip- and 
ip-. It is not yet clear what would be the probability in- 
terpretation for an action based on such Majorana-Weyl 
spinors. 

7. Antiparticles 

For every classical state r = [n+{x),n-{x)] we can define 
the "antistate" f = [n_|_(a;), n_ (x)] with n± = l — n±. The 
mapping r — > f transforms occupied bits into empty bits 
and vice versa. In the language of fermionic excitations it 
exchanges particles and holes or particles and antiparticles. 
A state with particle numbers (n+ , n_ ) is mapped into a 
state with — antiparticles of type + and n_ = n_ 
antiparticles of type — . For L space points Xi a state with 
n+ antiparticles has L~n^ particles of type +, and similar 
for n_ . For every state the number of particles plus holes 
for each species equals L, n+ + n+ — L , n_ + n_ = L, 
since for every space point there is either an empty or an 
occupied bit. 

The state with one antiparticle of type +, i.e. n_) = 
(1,0), can be written in the Grassmann formulation as 



91q = / Qi.oix)9lQ{x), 



(239) 

with gl Q — gi.o the Grassmann element conjugate to gifl, 
obeying eq. ([TUl) . 

I?#i,o(a;)5i,o(2/) = S{x - y). (240) 

One finds the simple expression 

gl^{x)=,k+{x). (241) 
Reordering the Hamiltonian (I195p . 

f ( d d 

n = -ij <^^d^ij+{x)-^{x) - d^ip-{x)-^{x) 

'(242) 

the wave function for one antiparticle obeys the same evo- 
lution law as for a particle. 



(243) 



Both the particle and the hole of type -I- are left-movers. 
Similarly, the wave function for one antiparticle of type 
— obeys 9t<7o,i — "dxQo^i such that both the particle and 
the antiparticle of type — are right-movers. 

More generally, a state with rl,+ antiparticles of type -I- 
and n_ antiparticles of type n_ can be written in analogy 
with eq. (|215l) as 



(n_|_!?T,_!) 



xi...XA,,yi...yfi 



with 



^,fi_(a;i ...Xfi^,yi...yfi_)^ Scgn+,n-{xi ■■■yn-) 

Sci'-iyn^)- . .-0-(2/i)V'+(a;ft+) ...i^+ixi), (245) 



1, 4 mod 4 and Sr 



-1 for 



whereby Sc = 1 for n+ -I- n 
•h^ + n_ =2,3 mod 4. For — n^, n_ — n_ the time 
evolution of qfi_^__n_ and qn+,n_ are given by identical (real) 
evolution equations. 



dtqn+,n.ixi-- 
dtqn+,n^ixi . . 



• J/n_ ) — Kji^ ,n_ qfi^ {xi . . . yn_ , 
■ yn- ) — -^n-i- ,n_ qri-^ ,n_ {xi . . . , 



(246) 



We define the operation of the charge conjugation C on 
the level of the Grassmann algebra g[ijj] (without ip) as 



Cg„+,n-{xi . . .y„_) = gn+,n-ix 



1 • ■ 



■Vn 



= 1, (247) 



implementing a mapping between particle and associated 
antiparticle states. The Hamiltonian commutes with C, 



(248) 



where we observe that for any Grassmann operator 



dtp 



the charge conjugate operator A'^ ~ 

CAC obtains by exchanging a and or -0 and d/di}} (keep- 
ing the order). For example, in a discrete formulation the 
particle number operators ()197p transform as 



Ml ^CN±C = L-N±. 



(249) 



The charge conjugate of the vacuum state |0) (|20ip is the 
unit element |1) — 1 of the Grassmann algebra, C|0) = 
|1),C|1) = |0), and obtains as 

qI+mAxi ■.■yn-)= Cgn^^n_{xi ...yn-) 



c 



d 



zC.C- 



d 



dtp+{xi) ■'■ dtj)-{yn_) 
■0+(a;i) . ..il)^{yn_). 



CC|G) 



(250) 



On the level of the wave functions the charge conjugation 
maps qT- q^ = {Cq)T and we can formally write 



(251) 



8. Complex structure 

The evolution equation (j246l) consists of two real equations 
for the real functions q and g'^, which have the same argu- 
ments. We can combine them into one complex equation. 
For this purpose we define the complex wave function 



_ 1 r 



Q-n+,n- +«(?«+,«_ -qn+.n-)]- 

(252) 

The real part of c is even under charge conjugation, while 
the imaginary part is odd. In this complex basis we can 
represent the charge conjugation as complex conjugation 



CCrt 



(253) 



The evolution equation (|246p becomes now a complex equa- 
tion 

dtCn+,n-{xi ...yn-)= KCn_^.n_{xi . ..yn.). (254) 
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The normalization of c„ 



is chosen such that 



XII. CONCLUSIONS 



(255) 



More generahy, a complex structure within a Grassmann 
algebra is defined by an involution 0,0^ = 1, which maps an 
arbitrary element of the Grassmann algebra g into another 
one, d{g). In our case the involution is defined by the 
charge conjugation, 0{g) ^ g" ^J^r (IrQr = J2t (1t9t- The 
presence of allows a map from a real Grassmann algebra 
to a complex Grassmann algebra, Qt ^ c-r- The elements 
which are even with respect to are mapped to Re{cr), 
while 0-odd elements are mapped to the imaginary part 
Im{cr), 

MCr) = \ilr+ Qr) , ^"^(Cr) = ^(^r - 9^)- (256) 

For the complex representation of the Grassmann al- 
gebra we only need one half of the basis elements gr- 
Indeed, from the knowledge of the complex coefficients 
Cn+,n_ we can extract both the real coefficients qn^,n_ and 
— ±'?L-n+,-L--n_ • We therefore need only the basis 
elements gn+,n^ with + n_ < L, since the associated 
Cn_|_.n„ also specifies qn^,n_ with L < n+ + n- < 2L. For 
n+ -|- ri„ — L we retain as independent basis elements the 
ones with n+ > L/2. For the special case n+ = — L/2 
- for even L - the charge conjugation flips between empty 
and filled positions and we only need half of them as in- 
dependent elements. With this restriction to independent 
Grassmann elements the complex Grassmann algebra (with 
complex wave function Cr) and the real Grassmann alge- 
bra (with real wave function qr) are equivalent. Each real 
Grassmann algebra admits the charge conjugation as a nat- 
ural possible complex structure. 

In the complex formulation we can exploit the structure 
of multiplication with complex numbers or the complex 
Fourier transform. In particular, the multiplication with i 
corresponds in the real formulation to a map qr q'^^q'^ ^ 
-Qt, 



1 

-1 



-1. 



(257) 



This is not equivalent to the formal multiplication with i in 
the real formulation of the Grassmann algebra. The latter 
has been employed above in order to illustrate the anal- 
ogy to the quantum formalism even though all equations 
remain real equations. (Only the real operator K = —iH 
appears in the evolution equations.) In contrast, whenever 
we work with a complex Grassmann algebra the multipli- 
cation with i is meant in the sense of eq. (|257l) . Using 
the complex formulation we recover the standard complex 
wave functions for multi-fermion states. 



We have established the map from a Grassmann func- 
tional Glipjip] to a family of probability distributions of a 
classical statistical ensemble {prit)} at different times t. 
In particular, this map also specifies the unitary time evo- 
lution of {prit)}. The existence of this map is very general 
- it only requires a relation between the "boundary terms" 
gf and gin in eq. (|M|) and a suitable form of the action 
S. Our examples cover two-state quantum systems as well 
as quantum field theories for fermions. "Diagonal observ- 
ables" find a direct interpretation in terms of the classi- 
cal probabilities. Such diagonal observables are realized as 
diagonal quantum operators acting on the real "classical 
wave function" q-r — Sr-\/Pr, = ±1, according to 

iAq)r = ^rplp^ ^rp = Ar6rp. (258) 

p 

They obey the standard classical rule for expectation values 

(259) 

r 

where At is associated to the value of A in the state r. Eq. 
(|259p is equivalent to the "quantum rule" 



(260) 



{A) = {qAq) = qrArpqp 

T T 

We also can express (A) as a standard functional integral 



(A) = / 2?VI?V^[V',V]G[V',V], 



(261) 



with A[ij},ip] a Grassmann element associated to the ob- 
servable. 

In general, the map G — >■ {prit)} is not invertible. While 
we can establish an invertible map between G and the fam- 
ily of real wave functions {q-r{t)}, we observe that two wave 
functions differing only by signs Sr will yield the same prob- 
ability distribution {^^-(t)}. This "gauge ambiguity" of the 
choice of the sign function Sr{t) does not affect the expec- 
tation values of diagonal observables (|260p . Furthermore, 
the choice of Srit) can be largely fixed by requirements of 
continuity (and differentiability) of {qrit)} in time (and 
space if appropriate). The residual gauge freedom can be 
fixed by a suitable convention which does not infiuence the 
prediction of expectation values. After gauge fixing, we 
can construct the inverse map {prit)} — >■ {qrit)}, which 
allows us to compute the wave function {qr} for a given 
probability density {pr}- For a suitable time evolution of 
{Pr{t)} we can map our classical statistical ensemble to the 
quantum formalism with wave function {qr(t)}. In turn, 
this can be mapped to a Grassmann functional G[ip,ip] and 
the associated functional integral. In this way we can ex- 
plicitly construct a quantum field theory of fermions from 
a classical statistical ensemble. 
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If the time evolution of the real wave function {qrit)} 
is compatible with a suitable involution defining a com- 
plex structure it can be mapped to the usual complex wave 
function of quantum mechanics. We have discussed several 
examples for such complex structures. In particular, the 
quantum field theory for two-dimensional fcrmions admits 
a complex structure associated to the notion of antiparti- 
cles, with involution corresponding to charge conjugation. 
We emphasize that different complex structures, based on 
different involutions are possible as well. For example, com- 
plex structures can be associated to discrete symmetries 
arising for quantum field theories with different "fiavors" 
of fermionic fields. 

The construction of a quantum field theory for fermions 
from a classical probability distribution is not limited to 
two dimensions. For example, we may generalize the ac- 
tion (I220p to four dimensions, where a real representation 
of the four matrices 7^ as 4 x 4 matrices exists. This in- 
volves instead of (x) four species of Grassmann variables 
ip.y{x). A real one particle wave function describes a Ma- 
jorana spinor with four components. In four dimensions, 
Majorana spinors are equivalent to Weyl spinors [l^. Ex- 
tension to a complex Grassmann algebra similar to sect. 
IXII accounts for Dirac spinors. One further may consider 
several flavors of fermions and add mass terms and inter- 
actions. 

Once a quantum field theory for fcrmions is obtained 
from a classical statistical ensemble, we may consider quan- 
tum states of a single particle or of several particles. Par- 
ticles can be interpreted here as isolated excitations of a 
vacuum state. Their properties depend on the choice of 
the vacuum state. We emphasize in this context that the 
state |0) in eq. (|20ip is not the unique possibility for a 
vacuum state. One could equally well consider the charge 
conjugate state C|0) or a linear combination of |0) and C|0). 
Another interesting possibility is the half filled state with 
an equal number of particles and antiparticles or holes. 

For the description of a particle state much less informa- 
tion is needed as compared to the whole classical prob- 
ability distribution {pr}- This is rather obvious, since 
{Pt} contains additional information for all multi-particle 
states. An isolated subsystem as a one particle state is 



therefore characterized by a coarse graining of the infor- 
mation. For our two dimensional example with vacuum |0) 
all one-particle observables must be computable from the 
information contained in the wave function qi^ (x) (|200l) or, 
more generally, from the associated density matrix. This 
generalizes to arbitrary isolated one particle excitations. 

In presence of a complex structure the complex one- 
particle wave function 'iPq\x) in quantum mechanics 
is typically constructed from two real "classical" wave 
functions qi{x) and q2{x),'ipQ\x) — qi{x) + iq2{x), cf. 
eq. (j252[) . The quantum mechanical probability density 

PqI^) = IV'Q''(a^)P = qlix) + 92(2^) is iiot enough to de- 
scribe the physics of interference. In our case, however, we 
dispose of two separate "coarse grained" probability distri- 
butions pi{x) = qf(x) and P2ix) = q^(x). Together with 
sign-information for qi and 52 this constitutes sufficient 
information for the description of phases and interference 
phenomena. Derivative operators as the quantum mechan- 
ical momentum can be associated Q to suitable classical 
observables Ar of the "microscopic classical ensemble" with 
states T. 

Often the classical correlation functions for the one- 
particle-observables can no longer be computed from the 
one-particle density matrix. In this case we have to deal 
with "incomplete statistics" [l^. The correlations be- 
tween measurements of different one-particle-observables 
have to be described by new types of correlation functions 
These new types of correlation functions can be im- 
plemented within the classical statistical ensemble, based 
on product structures for observables that differ from the 
classical product. In this way non-commutativity arises 
from the coarse graining of the information [8] . The neces- 
sary use of a non-commuting product for the description of 
measurement correlations in isolated subsystems explains 
[l|, Q why those violate Bell's inequalities . In this way 
the "no go theorems" [isl . [l6j for an implementation of 
quantum mechanics within classical statistics are circum- 
vented |ll|. Our examples of fermionic quantum field theo- 
ries constitute indeed explicit realizations of the emergence 
of quantum physics from classical statistics. 
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